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Abstract 



For simple inflationary models, we provide a consistent and complete scheme by 
which the macro-physical details of early universe inflation may be determined 
explicitly from the underlying micro-physical theory. We examine inflationary 
dynamics within the context of a fully renormalized, non-perturbative, and non- 
equilibrium quantum field theory using the closed time path formalism. We study 
the non-perturbative dynamics using the self-consistent Hartree approximation 
and the large limit. In addition, we include gravitation dynamically by means 
of the semi-classical approximation to Einstein gravity, allowing the cosmological 
geometry to be determined self-consistently by the evolution of the particle fields. 
We develop a simple and numerically implementable regularization and renor- 
malization scheme for which the renormalization counterterms contain no explicit 
dependence on the initial state. 

We concern ourselves first with the dynamics of new infiation models, showing 
how the dynamics become quantum fiuctuation driven. We then describe how these 
quantum fiuctuations may be reassembled into an effective field which behaves 
classically and acts as the source for the gravitational background. This allows 
us to determine the spectrum of primordial density perturbations in a consistent 
manner. We also examine explosive particle production after inflation through 
the processes of spinodal decomposition and parametric amplification, showing 
that the non-perturbative backreaction on the infiaton field prevents catastrophic 
particle production. We find a strong dependence of particle production on the 
expansion rate and on the symmetry properties of the infiaton. 



Preface 



Cosmology, the study of the beginnings of the Universe, has been a subject of 
study and speculation for thousands of years. We now live in a remarkable era 
in which cosmology has become an analytical subject of research to be studied, 
observed, and understood by science. 

Observations of the universe today have led to a description of the universe back 
to times of less than a tiny fraction of a second after the beginning of time. The 
most important advance in the study of our universe has been the construction of 
the Standard Cosmological Model. In analogy to the Standard Model of Particle 
Physics, this Model provides the framework around which all known aspects of 
Cosmology fit. 

The Standard Cosmology states that we live in an expanding universe in which, 
on the largest of scales, each point is like any other. While the details of the 
cosmology, such as the exact size, shape, and total content of matter, are yet 
to be determined, the uniformity of our universe has been verified to remarkable 
accuracy by relatively recent observations. 

The remarkable success of this model leads one to ask a very simple question: 
Why is the universe so uniform? It is an incredible fact that science has arrived at 
an answer on the basis of what is known about high energy physics and gravitation. 
The answer: Inflation. 

The solution to such a profound question deserves scrutiny and attention. The 
purpose of this dissertation is to examine what is known about inflation, to question 
its assumptions, and to carefully analyze its consequences and predictions. 
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Chapter 1 
Introduction 



The ultimate goal of studies of the early universe is the abihty to accurately predict 
the macro-physical properties of the universe from a micro-physical description. 
This is a daunting task which requires the relation of the full dynamics of the 
underlying quantum field theory to the geometrical dynamics of gravity. This 
evolution must then be followed over an extraordinary dynamic range before the 
results may be interpreted to yield concrete predictions for observable properties 
of today's universe. 

Such a program requires a number of features: 

• The ability to evolve quantum fields dynamically without any assumptions 
of an equilibrium state. 

• The ability to treat states of extremely high (non-perturbative) energy den- 
sities. 

• The ability to evolve gravity in a manner self-consistent with the evolution 
of the quantum fields. 

• The ability to interpret the results of the quantum field theory dynamics. 
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In addition to these fundamental features, one must be concerned with the renor- 
mahzation aspects of the field theory and also have means to explicitly evaluate 
the dynamics through analytical and numerical means. 

For a class of simple models, we have developed such a program. In the following 
pages, we apply this program to inflationary cosmology. 

1.1 The Need for a Non-equilibrium Description 

Non-equilibrium field theory [|l|, which describes the dynamical evolution of quan- 
tum fields, has become an important tool in particle physics and cosmology. The 
ability to compute observables in highly dynamical situations in which the results 
of ordinary equilibrium field theory fail, is necessary in order to gain a better under- 
standing of a wide number of phenomena, including phase transitions in condensed 
matter systems and at particle accelerators 0. Perhaps the most important ap- 
plications of non-equilibrium field theory, however, are to phenomena of the early 
universeil, 

The reason it is necessary to go beyond traditional equilibrium techniques in 
describing cosmology is the extreme environment of our universe at early times [§, 
^. It is now well accepted that early in its history, our universe was very hot 
with energy densities approaching the Planck scale, while the universe was rapidly 
cooling due to expansion. As a result of the rapid cooling, it is expected that 
matter goes through a number of phase transitions before we reach the present 



state of our universe[|T^, [TT], 0, |T^. 

Unless the matter fields involved in these phase transitions interact at a rate 
much greater than the cooling rate of the universe, the evolution proceeds out 



of equilibrium. This non-equilibrium evolution may in turn have a number of 
important physical consequences. One example is baryogenesis ||14|| . If one assumes 
that the universe began in a state with equal numbers of baryons and anti-baryons, 
then it is necessary to have some process in the history of our universe which results 
in a biased production of baryons over anti-baryons. Such a process, which may 
have occurred during the electroweak phase transition or at an earlier time, requires 
a departure from equilibrium since baryons and anti-baryons are produced in equal 
numbers in any equilibrium process. 

Another example of a non-equilibrium process in the early universe is the pro- 
duction of topological defects during a phase transition Just as the domain 
structure in a ferromagnet is very different when the sample is cooled adiabati- 
cally through the Curie temperature compared to a sample which is rapidly cooled, 
the production of topological defects depends strongly on the cooling rate of the 
universe [|16|. Again, equilibrium analyses will not provide an accurate picture of 
defect formation. 

The phenomena to be studied here, inflation O, Il8|, |T9|, EO, ^^ and the subse- 



quent phase of parametric amplification p2|, [23| , p^ , are likewise phenomena of a 
non-equilibrium nature [^, ^7\, Inflation, and the near exponential 

expansion which characterizes it, can only occur if the state of the inflating uni- 
verse differs from its equilibrium state. Furthermore, it is the dynamical evolution 
which allows inflation to end. Parametric amplification similarly relies on the rapid 
evolution of a quantum field to produce a parametric resonance of field quanta. A 
field in equilibrium shows no such behavior. 

Despite the non-equilibrium nature of these phenomena, the bulk of the work 



treating the processes has utihzed equihbrium techniques, most notably the effec- 
tive potential |]31]]. The effective potential is the name given to the equilibrium free 



energy of a system and while it is a very useful quantity in determining a number 
of important properties, such as the equilibrium ground state of the system, it 



is ineffective in describing the dynamics of the system[^]. It therefore becomes 
necessary to move beyond equilibrium techniques and to analyze the dynamics of 
these phenomena using non-equilibrium field theory. 

1.2 Previous Research in Non-equilibrium Quan- 
tum Field Theory 

With the recognition of the importance of non-equilibrium techniques, a signifi- 
cant literature has developed which deals with techniques, approximations, and 
applications of non-equilibrium field theory as applied to the early universe. 

The work of Hu and Calzetta and their collaborators]^ has focused primarily 
on the formal issues of non-equilibrium field theory, with an emphasis on statis- 
tical field theories. The Langevin approach derived from the influence functional 
has been shown to be a useful tool in studying the early universe[||. It should 
be noted, however, that statistical field theory remains an approximation to the 
full quantum field theory. Other important contributions of this research group 
include the study of the non-perturbative expansion for models with overall 
0{N) symmetry. In the context of the closed time path treatment of full quantum 
field theory, this work has included an analysis which places the large and other 
approximations within the context of the computation of the full two-particle irre- 
ducible (2PI) Green's functions P^l, 53, 53. Much of their work has been completed 



in expanding spacetimes and has included the issues of renormahzation as well as 
important applications. Of particular interest here is the work of Ramsey and Hu 
which treats the 0{N) theory in Friedmann- Robertson- Walker (FRW) spacetime 
including renormahzation [0^ , and applies the formalism to post-inflationary para- 
metric resonance p5| , reproducing and verifying many of the results of Refs. p9(] , 
which are included here in Ch. 5. 

Another research group with important contributions to the subject is the Field 



Theory group at the Los Alamos National Laboratory[^, ^ Their con- 
tributions include work on the expansion, including a determination of the 
next-to-leading order contributions in Minkowski spacetime as well as work on 



the leading order approximation in FRW spacetime |^8|. Other important con- 
tributions are provided by their non-equilibrium study of symmetry breaking ||39||. 
Also of importance to the present work is their analysis of the non-perturbative 
Hartree and large N approximations as compared to exact results in a toy quan- 
tum mechanics problem, which provides some indication as to the reliability of the 
approximations we use in our analysis pO|. 



The research which forms the foundation of the present work, however is that 
of Boyanovsky, de Vega, and Holman and their collaborators 2S, 41]. This 



important work includes the non-equilibrium study of symmetry breaking] 41 1 and 



domain formation in scalar field theories] 28] and also includes important contribu- 



tions to the formalism of non-equilibrium quantum field theory and its renormal- 



ization in FRW spacetime] 26]. In addition, the ground breaking work in collab- 



oration with Lee and Singh on reheating ][42|] is an important piece of the puzzle 



which along with the present work, paints a complete picture of inflation from its 
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earliest stages to its end, and on to the subsequent phases of the universe. 

Of course, all of the recent work in this area is grounded in the early pioneering 
work on non-equilibrium field theory [|I| and the study of quantum fields in curved 
spacetime|^^ which will be reviewed in further detail in Ch. 2. 



1.3 Content of the Thesis 

In this work we develop the techniques of non-equilibrium quantum field theory in 
conformally flat FRW spacetime, including non-perturbative approximation tech- 
niques and the introduction of a novel renormalization scheme which provides a 
simple interpretation of the renormalizations and allows for numerical analysis 
of the dynamics of quantum fields in both fixed and dynamically evolving FRW 



spacetimes[^, ^ |30|. We then provide two new applications of these tech- 
niques to important phenomena of the early universe: new inflation p7|, |30| and 
post-inflationary parametric amplification p6|, |29| . 



We begin our study in Ch. 2 with a review of FRW spacetime [§, ^, Q and the 



basics of field theory in an expanding universe |^3|. We describe the observational 
support for an FRW description including evidence for the important properties 
of homogeneity and isotropy and also provide the basics of the geometry of such a 
universe as described by general relativity. Next, we examine potential shortcom- 
ings of this standard picture, in the horizon, flatness, and relic problems, before 
moving on to the concept of inflation. Here, we describe how inflation provides 
solutions to the problems of FRW cosmology and discuss the elements necessary to 
a successful implementation of inflation. As we shall see, this will lead to the intro- 
duction of quantum field theory which in turn provides the particularly interesting 
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prospect that the same elements yielding inflation may also result in the very met- 
ric perturbations from which the large scale structure of our universe formed. We 
then introduce the basic concepts of field theory in FRW spacetime, including the 
introduction of conformal coordinates and a short discussion of renormalization in 
FRW spacetime. Finally, we complete the chapter with a short discussion of the 
formalism and physics of metric perturbations p5|, |46| . 

We follow our discussion of field theory in FRW spacetime with an in depth 
study of the techniques of non-equilibrium field theory in Ch. 3. First, we introduce 
the basic formalism before detailing the models and approximations we will find 
useful in our studies of the early universe. We pay particular attention to the non- 



perturbative Hartree|^, ^ and large A^[^, ^ 0, |3^, 0, ^ approximations, 
deriving the appropriate equations of motion and describing in detail the important 
properties and limitations of these tools. Next, we examine the initial conditions 
which will be used to begin the evolution of the dynamics when we apply our 
formalism to practical problems. This oft neglected detail will be seen to be of 
particular significance to both formal and practical issues of non-equilibrium field 
theory in FRW spacetime [^. The renormalization of the field theory is another 



issue of importance, especially when it will become necessary to numerically solve 
the equations of motion in a particular model. We show in detail how we arrive 
at our renormalization scheme and discuss the usefulness of our approach, while 
pointing out the important connection between the initial conditions for the field 
theory and its renormalization. Finally, we provide a description of the semi- 



classical approximation to gravity] 43], which allows one to determine the dynamics 



of the scale factor describing the classical expansion of the universe self-consistently 
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in terms of the evolution of the quantum fields of that universe. Such a description 
will also require us to make a number of additional renormalizations. 

Having dispensed with the details of the formalism of non-equilibrium field 
theory in FRW spacetime, Chs. 4 and 5 deal with applications of the techniques 



developed in Chs. 2 and 3. First, we examine the process of new inflation [|52 



53| , [27| , pO|] , beginning with a discussion of what the basic scenario is, why it was 
introduced, and what the problems with the scenario are. This leads into the 
introduction of the particular model of new inflation we study as well as a discussion 
of initial conditions. Next, we provide the analytical results of the dynamics of 
new inflation, followed by the details and results of our numerical study p7|, pO[| . 



The chapter is completed by a discussion of the dynamical assembly of a classical 
object which may be interpreted as driving the inflationary phase, followed by the 
details of the formation of metric perturbations and the de-coherence properties 
of the quantum field theory ||30|1. 

Our second application, discussed in Ch. 5, is the analysis of the phase of 
explosive particle production which may have occurred after a phase of inflation 
23, g, ER pi m m mm,m,mmm,m. We begin with a description of the 



basic scenario of chaotic inflation. We then analyze within our models the dynamics 
of phases of particle production due to the phenomena of spinodal decomposition 
and parametric amplification, including both analytical and numerical results [p^. 



29[| . We conclude with a discussion of the significance of our results in the context 
of reheating. 

Finally, we present our primary conclusions in Ch. 6, discussing the importance 
of this work to our understanding of cosmology and its interactions with particle 
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physics. We also discuss a number of other apphcations such as baryogenesis and 
defect formation which may benefit from a fuUy non-equihbrium analysis. 
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Chapter 2 

Field Theory in FRW Spacetime 



2.1 Introduction to FRW Spacetime 

The so-called standard cosmology presents a remarkably successful picture of our 
universe P, It is based upon two primary assumptions: 

1. The universe is homogeneous and isotropic on large scales. 

2. The geometry of the universe is determined by general relativity. 

Taken together, these properties yield a Friedmann-Robertson- Walker (FRW) cos- 
mology @ 0, 



The conditions of homogeneity and isotropy are very restrictive. If we define 
hij to be the metric on the three-dimensional surface defined by the spatial sections 
on the full four- dimensional spacetime, such that the differential spatial distances 
are given by 

dP' = hijdx'dx^ , (2.1) 



"'^This definition of an FRW cosmology is somewhat more general than that originally described 
by Friedmann [|63| , for which he made the additional assumption that the matter content of the 
universe is a pressureless, ideal fluid. While the properties of homogeneity and isotropy require 
the matter to act as an ideal fluid, we will allow for a more general equation of state. 
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then we can define the spatial Riemann curvature tensor R^fj^^. Given the connec- 
tion corresponding to the metric hij 

the Riemann tensor is 

The property of isotropy, which in fact, includes the property of homogeneity if 
required at every spatial point, requires that the spatial curvature defined in this 
way is constant. That is. 



'■J 



S-S- - S-^i] , (2.3) 



where the inverse metric k^^ has been used to raise indices. 

The constant /C may take on any value, but there are just three different cases 
corresponding to different topologies for the universe. If /C is negative, then the 
spatial sections take on a hyperbolic form. These may be defined as the three- 
dimensional surfaces embedded in a four dimensional Minkowski spacetime given 
by the relation 

J.2 2 2 2 7-)2 

t — X — y — z — R , 

where R is an arbitrary constant. The second case is /C = 0, which means that at 
every point in the space, the curvature is zero. The spatial sections are fiat and 
are therefore given by the three dimensional Euclidean metric. Finally, K, may be 
positive, corresponding to spatial sections of spherical form. We may define these 
by embedding them in a four dimensional Euclidean space. They are given by the 
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relation 



Again, R is an arbitrary constant. 

These three cases may be conveniently combined in the full four dimensional 
spacetime metric 

ds^ = g^^dx^dx" = dt^ - a^{t)dl^ , (2.4) 

— * 

where dl is the differential element giving the spatial sections and the scale factor 
a{t) is, for now, an arbitrary function of time. The spatial metric may be written 
in spherical coordinates as 

rfp = + r^dO^ + sm\e)d(P^ . (2.5) 

1 — Kr^^ 

Here, K = /C/|/C| takes on the discrete values K = —1, 0, +1 corresponding to the 
cases of hyperbolic, fiat, and spherical spatial sections respectively. We see that in 
the flat case of K = 0, the spatial metric may be written in Cartesian coordinates 
as 



The metric of the form given by ( p.4[) and ( p.5[ ) describing a homogeneous 



and isotropic cosmology is referred to as a Robertson- Walker metric] 64, 35 1. For 



a formal derivation of this important piece of the standard model, we refer the 



reader to the book by Wald||4^. Further discussions are given in the books by 
Weinberg ||6^ and Peebles 0. The latter reference also provides a nice discussion of 
the many equivalent coordinatizations of the Robertson- Walker metric including 
conformal coordinates, which are of particular importance and will be discussed in 
Sec. 2.4. 
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The second of the primary assumptions of the standard cosmology determines 
the dynamics of the scale factor a{t) in terms of the matter content of the universe. 
General relativity p^, ^ is based upon Einstein's equation relating the spacetime 



curvature to the energy- momentum T^^. The full Riemann curvature tensor, -R^^^, 
for the four dimensional spacetime is defined in terms of the spacetime metric g^i, in 
complete analogy to the spatial curvature (|2.2|) . The connection on the spacetime 
manifold is defined as 

and the curvature tensor is 

per fj "per a -per , pa per po per 

^J.l'p ~ '-'i^'' fip '^M-'- !/p ~r fip^ au up^ ap ■ 

Appearing in Einstein's equation is the Einstein curvature tensor G^j^, constructed 
from the Ricci tensor 

Raf3 = Rapl3 

and the Ricci scalar 

With these definitions, Einstein's equation reads 

G^,u = Rpu - ^Qpu'R' = -SuGnTi,^ , (2.6) 

where the proportionality constant is Newton's gravitational constant, which 
in natural units is related to the Planck mass by the equation G^ = l/mp^. 

Given the metric ( p.4| ) with (p.5|), Einstein's equation reduces to just two inde- 



pendent equations of motion for the scale factor: 



^2 

Gn = —3— — 3— = —SttGnTc. , (2.7) 



1^ 



a a? 0? 



Note that the non-diagonal elements of both G^i, and T^,y are identically zero while 
the remaining diagonal elements simply reproduce Eq. ( |2.8| ). We therefore see that 
T^j^ must have the form of an ideal fluid with energy density p = Tq and pressure 
p = —Tl in order to be consistent with the homogeneity and isotropy conditions. 



We can make a number of conclusions from the equations of motion ( |2.7| ) and 
( p.8[ ). First, if we assume that p > 0, we see that a^/a^ > and the universe is 
necessarily expanding (or contracting)^. Second, if p > then d/a < and the 
universe decelerates with time|^ 

We should note that for ordinary matter and radiation, the energy and pressure 
densities, which satisfy the respective equations of state p = and p = p/3, are 
functions of the scale factor and therefore also of time. In particular, ordinary 
matter has an energy density inversely proportional to the physical volume and 
therefore changes with the scale factor as 

Pm{t) OC ^ . 

Similarly the energy density of radiation changes as 

1 

Pr{t) OC 



which may be interpreted as a change in photon (or other massless particle) density 
proportional to l/a^(t) with an additional factor of l/a{t) due to the redshifting 

of the photon frequency. 

^If K — +1, however, there is the possibihty of a single moment in time satisfying SttGjv/ci/S = 
k/a^ for which a/a — 0. 

•^In fact, the universe decelerates for any equation of state satisfying p > —p/3. 
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Given the properties of a dynamic, decelerating universe, if the universe is 
presently expanding, as we will provide support for below, then physical length 
scales were much smaller at earlier times, with the expansion rate increasing as 
one goes back. If one extrapolates such a scenario of expansion far enough back, 
then one arises at the conclusion that there is an instant of time at which all 
physical length scales are reduced to zero. This is the source of the concept of the 
big bang, which should not be interpreted as an explosion of matter into space as 
its name implies, rather it represents the rapid expansion of space itself from an 
initial state approaching zero size and infinite energy density. 

Another property is that if the geometry of the universe is flat or hyperbolic 
with K — or K — —1 respectively, then an expanding universe expands for- 
ever. However, in the spherical case with K = +1, then at some instant in time 
SttGnp/S — K/a? and the expansion halts after which the universe begins to con- 
tract. In this case, the universe eventually will return to a state approaching zero 
size and infinite energy density called the big crunch. 

A final property of the expanding universe described by the standard cosmology 
that will be of particular interest is the existence of particle horizons. These exist 
because particles can travel at most at the speed of light. In particular, the furthest 
physical distance a light signal may have traveled at any time is given by 



This quantity is called the horizon distance. As we shall see, the important point is 
that if the scale factor grows no faster than linearly in t as is the case in matter or 
radiation dominated cosmologies, then the horizon distance decreases faster than 




(2.9) 
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the scale factor as one goes to earlier times. This means that points which are in 
causal contact today may have been causally disconnected at earlier times. We 
will return to this important point in the next section. 



Observational Support for the Standard Cosmology 

The picture presented here of the standard cosmology has received substantial ob- 
servational support, which has elevated the simple model to the level of a standard 
model0. The support comes in three different areas: 

1. Support for the theory of general relativity, especially on large scales. 

2. Support for the expansion of the universe. 

3. Support for the large scale homogeneity and isotropy of the universe. 

We now review the most significant evidence supporting a standard cosmology 
description of our universe. Complete discussions may be found in Refs. |^, 



Somewhat surprisingly, the element of support for the standard model which 
is perhaps the least certain is the theory of general relativity, largely because of 
the difficulties of verifying the theory on large scales. Of course, the success of 
general relativity on small scales is well known, beginning with the prediction 
and verification of the precession of the planet Mercury []66|, |67|1 . To some extent. 



however, we must simply assume that general relativity is the correct theory of 
gravity on all cosmological scales. 



''^The terminology of the standard model of cosmology was introduced by S. Weinberglp^l in 
analogy to the incredibly successful standard model of particle physics. A standard model has 
come to mean a simple, yet complete, theoretical model which has held up under significant 
experimental or observational testing. 
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Observational evidence for the expansion of the universe, on the other hand, is 
very extensive. The expansion was first observed by Hubble ^8|, who found that 



light emanating from distant cosmological objects was Doppler shifted toward the 
red part of the visible spectrum. Such redshifting indicates movement of objects 
away from our galaxy. Furthermore, Hubble observed that the amount of redshift- 
ing increases with the distance of the object. In particular, Hubble determined 
that there was a near linear relation between an objects redshift 1 + z, defined as 
the ratio of the observed wavelength to the expected wavelength for that source, 
and the distance d. This relationship is expressed as Hubble's Law: 

where Hq is Hubble's constant. 

While deficiencies in the ability of Hubble to accurately measure the distances 
of cosmological objects led to a value of Hq an order of magnitude larger than recent 
determinations, Hubble's Law itself is now well verified^. We therefore know quite 
definitively that the universe is undergoing expansion. 

The final piece of evidence in support of the standard cosmological model is 
direct observation of the homogeneity and isotropy of the large scale universe. 
The primary and most convincing observation in this regard is provided by the 
Cosmic Background Radiation (CBR) which was found to have an almost perfect 
blackbody spectrum at a temperature T = 2.73K with little variation in this 
temperature with the direction of the observation. In fact, if one subtracts the 
dipole component of the temperature variation, which may be explained by the 



^This relation is in fact just the first term in a more complete expansion of the form z 
Hod + {qo — l){Hod)^ /2 + • • •. Hubble's Law is accurate for distances d < I/Hq. 
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movement of our reference frame with respect to the rest frame of the CBR, the 
temperature variation over the sky is of order 5T /T ~ 10"*^. This means that at 
the time the photons making up the CBR last interacted with ordinary matter, 
the universe was homogeneous and isotropic to one part in 10,000. To better 
understand what this means, it helps to understand the source of the CBR. 

Our standard picture of the early universe indicates that its temperature was 
once many orders of magnitude hotter than it is today. At high enough temper- 
atures, we know that electrons dissociate themselves from protons. At such high 
temperatures, the radiation in the universe continually interacts with the charged 
electrons and protons, and photons are in a state of thermal equilibrium. However, 
as the universe cools, the electrons combine with the protons during a time called 
recombination (despite the fact that these particles were never combined before 
this time). Without the charged particles on which to scatter, the photons become 
effectively non-interacting. As a result, they simply free stream until they reach 
observers today. 

In the meantime, the photons have redshifted with the expansion of the uni- 
verse, and since a redshifted thermal (blackbody) spectrum remains thermal, we 
observe this cosmic background radiation as a near perfect blackbody spectrum at 
a temperature given in terms of the photon temperature at recombination Trc 
by the relation 

'"(1 + -)- 

The redshift at recombination is approximately z ~ 1200. 
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Since the ratio 5T/T is a constant for the cosmic radiation background during 
its evolution from recombination to today, it is clear that the recent observations of 
an isotropic and homogeneous CBR indicate that the universe was very isotropic 
and homogeneous at the time of recombination. 

There are a number of other observations lending support to the model of an 
isotropic and homogeneous universe. The X-ray background, due to unresolved 
high redshift sources has been determined to be isotropic to about 5%. Similarly, 
the distribution of faint radio sources shows no preferential direction. Furthermore, 
the pecuhar velocity field of galaxies on scales 50-lOOMpc, which gives a direct 
measure of density variation through the movement of galaxies relative to the 
Hubble flow, provides a density contrast 5p/ p ~ 0.1. 

Finally, galaxy distribution itself is found to be generally isotropic, although 
it is somewhat disconcerting that we have observed structure at even the largest 
observable scales. If the universe is truly isotropic on the largest scales then at some 
point we should reach a limit beyond which there are no new, larger structures. 

2.2 Shortcomings of the Standard Cosmology 

As we have seen, the standard cosmology is very well founded. However, it is 
not without problems. Obviously, the universe is not perfectly homogeneous since 
there are a number of individual structures ranging from planets to superclusters. 
These presumably form over time through gravitational instability, acting on the 
initially small density variations we observe in the CBR. However, the standard 
cosmology does not account for these density fluctuations, and it is fair to ask 
about their origin. 
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One of the most interesting challenges to the standard cosmology, on the other 
hand, is the question of how the universe could be so homogeneous in the first place. 
This would seem like an empty question given our assumptions of homogeneity were 
it not for the concept of particle horizons. As we saw in the preceding section, the 
horizon size in a matter or radiation dominated cosmology decreases relative to 
physical length scales as one goes to earlier times. This means that the horizon 
size at the time of recombination was much smaller than it is today, and there 
would have been no way in which correlations in the GBR in different directions 
could have formed. 

In fact, each causal region at the time of recombination subtends less than 
a single degree of the presently observable universe and one might expect each 
of these causal regions to have evolved somewhat independently, resulting in a 
microwave background temperature spectrum which varies significantly from one 
region to the next. The actually observed homogeneity may, of course, be explained 
by simply assuming that the universe began in an almost exactly homogeneous 
state independent of particle interactions. However, we will see that there is a 
much more eloquent solution to the horizon problem based on a marriage of particle 
physics and cosmology. 

A third problem of the standard cosmology is the relative flatness of our current 
universe. From the Friedmann equation ( p.7|) , we see that there are two terms 
contributing to the expansion rate of the universe. The first is due to the energy 
density of matter and fields and the second is due to the curvature of the universe. 

25 



It is convenient to introduce the critical density: 

3 a2(t) 
^'^^^ = SttGm a' it) ■ 

If p = Pc then we necessarily have a critical, spatially flat geometry while p < Pc 
and p > Pc correspond respectively to the open, hyperbolic geometry {K = — 1) and 
the closed, spherical geometry {K = +1). This behavior is often summarized by the 
introduction of the single parameter Q(t) = p{t)/pc{t), with = 1 corresponding 
to the critical universe. 

While the present value is not well determined, through indirect measure- 
ment as a result of the successes of Big Bang Nucleosynthesis combined with a 
number of more direct measurement methods, it is restricted to the modest range 
0.01 < Qq < 2.0, while most cosmologists would agree that Qq is very likely to lie 
in the range 0.2 < l^o < iQ. 

If, in addition to Q, we define the quantity flR{t) = —K/a?(t), Friedmann's 
equation ( p. 7] ) becomes 

n(t) + nn{t) = 1 . (2.10) 

This equation displays a remarkable coincidence: given the accepted range of fl 
and assuming k ^ (that is, Q ^ 1), we see that the contributions of Q and 
flu to unity in Eq. ( p.lO| ) are within an order of magnitude or so of each other. 



What makes this so remarkable is that these two terms have a very different time 



^In fact the determination of fig is an exciting and fast moving area of observational research. 
The current status of this research is perhaps best iUustrated by two recent determinations. The 
first measures the acceleration of distant Type la supernovae, which are assumed to be standard 
candles, and favors values of fio < 0.5 while effectively ruling out the possibility of fio = 1 unless 
a cosmological constant is added p9||. However, another observation which uses number counts of 



distant x-ray clusters strongly favors fig > 0.5 70|. While progress is being made, there clearly 
is yet to be a consesus on fig- 
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dependence since oc Xjc? while oc l/[a^(i)a(i)] for a matter dominated 
cosmology and f2 oc l/[d^(t)a^(t)] for a radiation dominated cosmology. For these 
two quantities, which have drastically different time dependence, to be of the same 
order at this instant of time despite the huge variation of a{f) of many, many orders 
of magnitude during the history of our universe is extraordinary. This coincidence 
is referred to as the flatness problem since it does not exist if the geometry of the 
universe is flat since that would mean fi^j = 0. 

Another way of posing this problem is to note that since decreases faster 
with expansion of the universe than ^^^(t), the condition that these quantities now 
be within a couple of orders of magnitude of unity extrapolate back to the condition 
moments after the Big Bang of il(ipi ~ 10~^^s) ~ 1 ± 10"^'^. This indicates an 
extreme fine tuning of the initial state of the universe. 

The final problem of the standard cosmology arises as a result of the Grand 
Unified Theories (GUT's) of particle physics. These theories, which unify all the 
known, non-gravitational forces of the universe, typically predict the production of 
topological defects. In particular, all known GUT's predict the production of large 
numbers of monopoles which would quickly come to dominate the energy density 
of the universe. Since these objects are topologically stable, they would continue 
to dominate for all time, and the cosmology we observe today could not exist. A 
somewhat similar problem arises with the possible production of domain walls in 
early universe phase transitions, which would be expected to disrupt the isotropy 
of the universe. Together, such potential problems constitute the relic problem. 
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2.3 Introduction to Inflation 



It is a remarkable fact that through the marriage of particle physics and cosmol- 
ogy, a single mechanism has been discovered which addresses all of the problems 
discussed in the preceding section: Inflation. 

Inflation|r7| is a rapid expansion of the universe, characterized by acceleration 



of the scale factor: 

d{t) 



>0. (2.11) 



a{t) 

The proto- typical inflationary universe expands at an exponential rate with a/a = 
H with H constant. This is the so called de Sitter universe [[7I|]. 

We now describe how such a universe solves the horizon, flatness, and relic 
problems; we will reserve the fluctuation problem for later. We will use the ex- 
ponentially expanding de Sitter universe with a{t) = aoexp{Ht) as an example, 
although the arguments carry over, albeit to a lesser effectiveness, for any case 
satisfying ( p.ll| ). 



For the de Sitter universe, the horizon size (12.91) is a constant 



^ dt _ 1 



Clearly, is such a situation, the physical length scales decrease faster than the 
horizon distance as one goes back in time. For a concrete example, let us assume 
that the de Sitter phase lasts a time r. An initial region of radius r is therefore 
expanded to size r exp{Ht). If we take this region to be homogeneous and causally 
connected initially, then the expanded region will also be homogeneous. If we make 
r large enough, we can therefore inflate an initially small region of the universe 
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creating a homogeneous region which is large enough to fully contain our present 
universe. 

We therefore see how inflation solves the horizon problem. The process en- 
sures that distant regions of our observable universe which could not have possibly 
interacted during the matter and radiation dominated phases of expansion had in- 
teracted at a much earlier time before inflation took place. In this sense, inflation 
acts as a mechanism to set the initial conditions for the standard cosmology such 
that the present universe is homogeneous on large scales. Note that the horizon 
problem requires of the order of 60 e-folds of inflation to take place in order to 
produce a homogeneous universe of the present size. In the exponential case, we 
therefore require that Ht > 60. 

The simplest solution to the flatness problem is to remove the contribution of 
flu entirely from Eq. ( |2.10| ), such that = 1. Such a contribution may be removed 
by reducing the relative contribution of the curvature term in Friedmann's equation 
( p.7|) with respect to the energy density. The requirement to do so is to have a 
component of the energy density which decreases more slowly than l/a^(t). If 
this requirement is satisfied then sufficient expansion of the universe will reduce 
the contribution of relative to Q. Having the universe dominated by energy 
density satisfying this requirement is equivalent to the condition a/a > 0, precisely 
the condition for inflation. 

In simple terms, inflation solves the flatness problem by taking an initial region 
which is small enough to have negligible curvature and expanding it to the size of 
our present universe. We also see that in solving the flatness problem, inflation 
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provides a powerful prediction: Q = 

The relic problem is solved in a similar way to the flatness problem. As 
monopoles act as ordinary matter, their energy density decays with expansion 
at least as fast as l/a^{t). It is therefore clear that the same conditions required 
to solve the flatness problem will necessarily dilute the number of monopoles to a 
negligible level. 

If it is also necessary to dilute the energy density in domain walls, then this 
is also possible, although it may require a dominant contribution to the energy 
density which falls more slowly then the l/a(t) dependence of the energy density 
of infinite domain walls. 

Now that it is clear how inflation solves the horizon, flatness, and relic prob- 
lems of the standard cosmology, it is important to determine under what conditions 
inflation occurs. From Einstein's equations ( p.7|) and (|2.8|) , we see that the require- 
ment d/a > translates to the condition on the equation of state of the universe 
of 

p + 3p<0. (2.12) 

This requires that the pressure density of the universe is negative with p < — p/3. 

At first, such an equation of state (which breaks the strong energy condition||44||) 
seems odd. However, it is relatively easy to produce such a situation in field theory. 
As an example, let there be a scalar field in Minkowski space with Lagrangian 

L = d^<j)d^<j) - V{<j)) . 

^There are a number of open inflation models which aUow fl < 1 through a mechanism of 
false vacuum tunneling. However, these models are somewhat awkward in their construction and 
implementation. See [73| . 
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The energy density due to is 
while the pressure density is 

We see that if a situation arises in which the contributions of the potential V{4>) 
dominate over the kinetic {<p) and gradient (V0^) terms in these quantities, we 
arrive at an equation of state p = —p. Such is the case when is constant and 
different from its vacuum value. 

This example also provides us with an example of the importance of quantum 
field theory in producing inflation, since a hydrodynamic description of the universe 
containing classical matter does not allow an equation of state p < 0. 

It should be noted that the effect of a potential energy contribution to the 
energy density is similar to that of a cosmological constant A in Einstein's equation: 

Gij.u — ^g^u = —SiTtGnThi, . 

A cosmological constant dominating over the energy density of matter will there- 
fore produce inflation. However, even in this case, quantum field theory will be 
necessary to produce an end to inflation. This will be important in our study of 
new inflation in Ch. 4. 

Finally, we return to the one problem of the standard cosmology we have yet 
to address: the existence of small primordial density fluctuations. While inflation 
has the effect of homogenizing the universe, it surprisingly also provides density 
fluctuations on cosmological scales. 
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The basic mechanism, which we will describe in more detail in Sec. 2.5 after 
introducing the basics of field theory, is that inflation takes quantum fluctuations 
on microscopic scales and inflates them to cosmic scales. It is these quantum 
fluctuations which constitute the density fluctuations from which the structure of 
our universe formed. 



Quantum field theory in curved spacetime is a rich and interesting topic. Here we 
restrict the discussion to the areas of primary interest to the present study. We 
present the basic elements of scalar field quantization in FRW spacetime, providing 
the field equations and describing the mode decompositions of the field. We pay 
particular attention to the concept of conformal coordinates as this will be of 
importance in constructing our renormalization scheme. We complete this section 
with a short review of renormalization in curved spacetime. Extensions to fields 
of higher spin and to more general spacetimes are provided in the literature [Q. 
The Lagrangian for a general scalar field in curved spacetime is 



where is the usual covariant derivative and we allow for a potential V^($; x) 
which may depend explicitly on the spacetime coordinates due to the existence 
of curvature scalars such as the Ricci scalar. In FRW cosmologies, this explicit 
dependence is restricted to the time variable. We will be primarily interested in 
the theory, which represents the simplest, non-trivial, interacting scalar field 



2.4 Field Theory in FRW Spacetime 




(2.13) 
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theory. In this case, the potential has the form 



Vmx);t] 



1 r 



(2.14) 



Here, we have included a term which couples the scalar field directly to the scalar 
curvature 7?.(t) through the coupling constant ^. There are two values of ^ which 
receive particular attention: 

• minimal coupling with ^ = 0, and 

• conformal coupling with ^ = 1/6. 

The importance of the former is obvious from its simplicity, while the importance 
of the latter will be discussed when we introduce conformal coordinates. Of course, 
like any other coupling constants, without any strong physical principle to tell us 
the value of ^, we must assume that this constant may take on any value and allow 
it to be determined by measurement. 

The action for a scalar field in curved spacetime is simply 



where g = det g^„. Variation of this action with respect to the field $ provides the 
field equation 



where the prime denotes differentiation with respect to For the A$ theory in 
FRW spacetime, this becomes 




(2.15) 



(2.16) 



W'Vf. + m'^ + ^n{t) + -^'^{x) $(x) = 0. 



(2.17) 
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In order to seek solutions of these field equations, it is convenient to introduce 
a mode decomposition of the field Restricting ourselves to a free scalar field in 
an FRW cosmology, we may write 



dfi{k) 



(2.18) 



where and at are respectively the annihilation and creation operators for the 
field quanta. These operators satisfy the relations 





- ^kk' ' 


(2.19) 




= 0, 


(2.20) 


4' 4. 


= . 


(2.21) 



We may make a separation of the modes into their space and time dependent 
pieces 

M^) = yd^Wd^) ■ (2-22) 
The mode functions ^/^(t) satisfy the mode equations 

U^.it) = . (2.23) 



— + 3-7TV— + A; + + ^7^(t) 



dt'^ a{t) dt 
The spatial functions 3^^(a;) satisfy the relation 

V*V,3^.(f) = -{e - K)yAx 



k\-^J ' 



(2.24) 



where Vj is the covariant derivative on the three dimensional spatial metric hij of 
Eq. (^). 

The solutions to ( ^.24[ ) have been determined for each of the three FRW cos- 
mologies. In the fiat, K = case, they are simply 



yj^{x)=e"'-^ {K = 0) 
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(2.25) 



with k = {kx, ky, kz). For the closed and open cases respectively, these functions 
may be written as 

y^^{x) = IiUx)yf{e, <t>) (K = ±1) . (2.26) 

Here, k = (fc, J, M) and the Yf^ are the usual spherical harmonics. For K = +1, 
k takes on positive integer values with J = 0,1,..., A; — 1, while foT K = —1, k 
takes on any positive indefinite value while J is a positive indefinite integer. In 
either case, M takes on the values M = —J,—J + 1,...,J. 
The function n^j(x) is 

^uix) = [Ik'ik' + + J^)}"'^' sinh-^ X [j^) cos {kx) , 

(2.27) 

and n^j(x) is given by the same expression taking x —^X- 

Finally, we return to the measure dfl{k) appearing in the mode decomposition 
( ^181) . We have, 

/# (^ = 0) 

Em,m {K = +l) (2.28) 



j djl{k) 



JdkEjM (^ = -1) 

In what follows and throughout, we restrict ourselves to the case of a spatially 
fiat {K = 0) FRW cosmological model. This is done for convenience and because an 
initially curved space takes on the properties of a fiat universe soon after infiation 
begins. 

Conformal Coordinates 

Quantum field theory may be formulated using a number of different coordinate 
systems in which the field equations take on different forms, while the physics 
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described is the same. Of these, few have particular physical significance. However, 
the system of conformal coordinates in fiat FRW spacetime is of importance due 
to its relationship with Minkowski spacetime. 

A conformal transformation is one in which the spacetime metric is multiplied 
by an overall factor. That is g^^, —>■ C{x)g^,^. The new spacetime metric C{x)g^j^i, 
is said to be conformal to g^^,. If we examine the flat FRW metric 

ds^ = dt^ - a^{t) \dx^ + dy^ + dz 

we see that if we define the conformal time variable 

then the metric becomes 

ds^ = 0^(77) [drf" - dx^ - dy^ - dz^] , (2.29) 

which is conformal to the Minkowski metric. This is of particular importance 
because wc may define a vacuum in conformal coordinates in which the adiabatic 
limit d/a is equivalent to the Minkowski vacuum. This presents us with a 
foundation on which we can build our vacuum state. 

In addition to working in conformal time, it is convenient to make a conformal 
rescaling of the field by writing 

$(f,i)^^^i^. (2.30) 

The action in terms of the new field xi.^i'H) becomes, after integration by parts 
and dropping the surface term 

j d'xdri l^x' - I (Vx)' - V(x)} , (2.31) 
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where primes denote differentiation with respect to conformal time t] and 

with the Ricci scalar TZ = Qa' {rj) / a?{7]). Here we see that the action ( p.31| ) takes 
on the same form as that of a scalar field in Minkowski spacetime with a potential 

v(x). 

If ^(x, t) is a field with A^^ potential ( |2.14| ), then the conformal potential is 



v(x) 



2 I A 4 



2 

From this, we recognize an important property of the action ( p.31 ). If t) is 



a massless {m? = 0), conformally coupled = 1/6) field, then the action for the 
conformal field xi^^v) ( f^-31D is invariant under conformal transformations. 



Renormalization in Curved Spacetime 



Renormalization is a delicate issue in curved spacetime for a number of reasons]^ . 
First, whereas in Minkowski spacetime the zero point energy has no physical mean- 
ing, such an energy becomes important in the gravitational field equations. A di- 
vergent zero point energy therefore cannot be removed in a more or less arbitrary 
manner as in Minkowski spacetime, rather it must be subtracted in such a way to 
yield a sensible finite result. Second, one must be concerned about preserving the 
property of covariant conservation of the energy momentum tensor. 

We will require the following properties of our renormalization scheme: 

1. It must result in numerically amenable renormalized equations of motion. 

2. It must preserve covariant conservation of the energy momentum in the limit 
that the regulator is removed. 
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3. All subtractions must be removed through time-m(iej>en(ientrenormalizations 
of the parameters of the theory. 

4. The zero point energy must vanish in the true vacuum of the theory. 

The standard methods of dimensional regularization and covariant point split- 
ting are inherently covariant schemes. However, they are not well suited to nu- 
merical analysis. An exception to this rule is provided by the recently constructed 



renormalization scheme of Baacke, Heitmann, and Patzold[^ which uses dimen- 
sional regularization combined with a clever separation of the divergences from the 
finite terms to provide a numerically amenable and inherently covariant regular- 
ization scheme. 

We choose to turn to a non-covariant scheme for which we have to verify co- 
variant conservation after the fact. We introduce a large momentum cutoff to 
regularize divergent integrals, and show that the resulting cutoff dependent terms 
may be used to renormalize the parameters of the theory. The scheme satisfies all 
of the required properties and in particular allows for efficient numerical analysis 
of the resulting equations of motion. 

2.5 Metric Perturbations 

One of the most important aspects of inflationary theory is its ability to produce 
primordial metric perturbations on scales relevant to the formation of large scale 
structure and to the temperature inhomogeneities of the GBR. Here, we review the 
elements of the formalism of cosmological perturbations that will be important for 
our study of inflation. We utilize the gauge invariant formalism of Mukhanov, 
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Feldman, and Brandenberger|^5|, to which we refer the reader for a more complete 
and detailed exposition of the subject. 

Obviously the treatment of cosmological perturbations requires a departure 
from the purely homogeneous FRW cosmology. We therefore write the metric in 
terms of a background FRW spacetime denoted by ^^^g^u and a perturbation metric 
SQiiu which will be considered to be small: 

Of course, there are 10 independent components making up Sg^,^. Those may be 
categorized in terms of their transformation properties under spatial coordinate 
transformations. In particular, the perturbations may transform as spatial scalars, 
vectors, or tensors. In conformal time, we may writeQ 

2 / -20 Bii-Si \ 

Sg,. = -a (r^) . _ _2 _ j + . + p^^^ + j..^ j , (2-32) 

where hij is the spatial background metric and /|j is the covariant derivative 
with respect to hij. The vector components satisfy the constraints 



while the tensor component satisfies 



to yield a total of 10 independent functions. 

^Here, our notation differs somewhat from since we continue to represent the spatial 
metric by hij {'^ij in |E5[). We use fey for the tensor perturbation {hij in p5[). 
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There exists a gauge freedom in general relativity which corresponds to a choice 
of coordinates on the curved spacetime, and it is important to note that the met- 
ric perturbation is not invariant under these gauge transformations. A general, 
infinitesimal coordinate transformation is given by 

^ 5" = + r • (2.33) 

The transverse piece of the 3- vector ^* satisfying Cj^ii ~ contributes only to vector 
perturbations. We will be interested in the remaining pieces and ^ where the 
latter quantity is defined by the expression 

Given a gauge transformation which preserves the scalar nature of perturbations, 
with = 0, the four scalar functions of 6gf^u transform as 

4> = cf^-^e-e', (2.34) 

a 

iP = + (2.35) 
a 

B = E + (2.36) 
E = E-^. (2.37) 

From these, it is possible to construct two gauge invariant combinations. The 
simplest are 

1 r. 



$ = 0+ - 
a 



{B-E)a , (2.38) 



■qj = ^-^m-E') . (2.39) 
a 

There exists, of course, a gauge freedom corresponding to the transformation 
( |2.33 ). Of particular note is the longitudinal gauge which is defined hj B = E = 0. 
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As a result, 0/ = $ and ipi = where the subscript denotes that the quantity is 
evaluated in the longitudinal gauge. 

We will be interested in metric perturbations due to fluctuations in a scalar 
field (f{x,t), which we separate into a homogeneous, classical part v^o(''^) and an 
inhomogeneous fluctuation 6(p{x,t). We write 

where it is assumed that 6ip <^ (fo- For metric perturbations which result from 
the fluctuations of a scalar fleld, it turns out that ^ = $ and we may reduce 
the perturbation equations to include just one of these gauge invariant quantities. 
We therefore focus on the dynamics of $ which is referred to as the Bardeen 
potential]^. In the longitudinal gauge, $ may be identified with the Newtonian 



potential, the gradient of which provides the gravitational acceleration of a test 
particle relative to the global inertial coordinates. The importance of this quantity 
is due to its close relationship with the power spectrum per logarithmic k interval 
which gives the amplitude of scalar metric perturbations. For modes which 
have physical wavelength larger than the horizon scale, k < a{t)dh{t), this relation 
is simply 

where $fc are the Fourier components of $. 

The power spectrum |5fe(t)p is directly related to the temperature variation 
of the cosmic microwave background radiation on large angular scales [^] (greater 
than one degree). Assuming that the last scattering surface at recombination is 
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infinitely thin, this relation is simply 

AT{e) 

T 



(2.40) 



where 6 is the angular separation corresponding to the wavenumber kg at recom- 
bination. 



The equation of motion for these Fourier modes of the Bardeen potential isPS 



= . 

(2.41) 



Here, H = a/a is the Hubble expansion rate. The Bardeen potential may also be 
directly related to the field fluctuations 6ip, or rather the gauge invariant combi- 



nation 5^9^'* = 6(p + (po{B — E ), by the constraint equation |^ 

Air 



(2.42) 



This gauge invariant scalar field fluctuation obeys the equation of motion 

Vfc' - 40<i>fc + 2V'M<l>k = , (2.43) 

where A4 is the (generally time dependent) mass and V{ipo) is the potential for 
the classical field (p. 

While the above equations are useful in determining the power spectrum for 
modes of physical wavelengths greater than the horizon size, it is not until these 
fluctuations re-enter the horizon long after inflation has ended that they will be 
relevant to structure formation or observations of the CBR. A convenient way of 
progressing from the spectrum of fluctuations during inflation to the spectrum of 
density perturbations at re-entry is the conservation law|?3| 

, (2.44) 



2 ^k/H + 



3 l+p/p 
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which is vahd after the first horizon crossing of modes with wavevector k < adh- 



43 



Chapter 3 



Non-equilibrium Field Theory 

3.1 Introduction to the Closed Time Path For- 
mahsm 

The out of equilibrium dynamics of a quantum field is obtained from the Closed 
Time Path formalism of Keldysh and Schwinger[^]. The starting point for deter- 
mining the time evolution of a quantum system is the Liouville equation. This 
equation determines the evolution of the quantum density matrix, p[(/), in 
terms of the Hamiltonian operator H{t) of the system. In the Schrodinger picture, 
in which operators are constant while the states change with time, the Liouville 
equation reads 

'^ = m),m]- (3.1) 

Here, we allow for an explicitly time dependent Hamiltonian H{t), as will be nec- 
essary to treat quantum fields in the time dependent background of an expanding 
universe. 

A formal solution of this equation may be written in terms of the time evolution 
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operator 



4 J^^ dtH{t 



U{t,t) = exp 
in the form 

pit) = Uit,to)pito)U-\t,to). (3.2) 

The quantity p(to) determines the initial state of the system. One possible choice is 
to let the initial state be one of local thermodynamic equilibrium at a temperature 
T = The unnormalized initial density matrix then takes the form 

/3(to) = exp {-fSHito)) . (3.3) 

While it is not obvious that this choice is the most appropriate initial state for 

the highly non-equilibrium processes of early universe, it is a very convenient and 

simple choice which includes both the case of an initial finite temperature as well as 

that of an initial state at zero temperature (the vacuum state) in the limit (3 oo. 

The latter state may be particularly appropriate for the dynamics after an initial 

stage of inflation. Of course, there exist other possible initial states. 

Given the evolution of the density matrix ( p.2| ), ensemble averages of operators 

are given by the expression]^ (again in the Schrodinger picture) 

^ Tr[pit)0] ^ Tr[pito)Uito,t')Uit',t)OUit,to)] ^ 
Trpit) Trpito) 

where we have inserted the identity, U{t, t')U{t' ,t) with t' an arbitrary time which 

may be taken to infinity. The form of the initial density matrix ( p.3| ) allows us to 

extend this expression into the imaginary time plane by writing 

_ Trpito - tp, to)U{to, t')Uit',t)OU{t, to)] 
TrU{to-i/3,to) 



^From here on, we drop the hats from the operators. 
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Figure 3.1: The Closed Time Path contour in complex time. 



We give the following interpretation to this expression, reading from right to left. 
The state is first evolved forward from the initial time to to t when the operator O 
is inserted. We then evolve this state forward to time t', back again to the initial 
time, and down the imaginary axis to to ~ This is normalized by the evolution 
in the imaginary time direction from tg to to — il3. This process is depicted by the 



contour in Fig. This is the Closed Time Path (CTP) formalism. Notice that 
the expectation values defined in this way are in-in expectation values as both the 
initial and final (real) times appearing in the trace (|3.5| ) are given by to, and all 
observable quantities are determined with respect to the Hamiltonian at the initial 
time. 
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The actual evolution of various quantities in the theory are evaluated by ei- 
ther constructing the appropriate Green functions [^], or by choosing an explicit 
ansatz for the functional form of the time dependent density matrix so that the 
trace in (|3.4|) may be explicitly evaluated as a functional integral[^]. Within the 
approximations we will use, the methods are equivalent. For completeness, we will 
describe both methods, beginning with the Green function approach. We focus 
on the case of a real scalar field with no derivative couplings appearing in the La- 
grangian. The analysis will be generalized to the case of a vector 0{N) multiplet 
in a later section. 



Computation of the Green Functions 

We introduce the generating functional with sources, J"*", J~, and J^, defined on 
each of the three legs of the contour in Fig. lO 



Z[J+, J~, J^] = Tr \u{t, - to; J^)f/(to, t ; J-)U{t , to; J' 



(3.6) 



where we have defined the sourced time evolution operators which for a scalar field 
$ have the form 



U{t, t';J) = exp |-i dt" !^H{t") - J d^xJ{x)<l>{ 



X 



(3.7) 



By inserting a complete set of field eigenstates, which we shall denote as $, between 
the time evolution operators, this may be written as 



Z[J+, J-, Jf^] = J V<!?+V<^-V<!>f^ exp |z J d^x ^* dt{£[$+, J+] - J-]} 

(3.8) 



to-i/3 



Ito 
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with the boundary conditions $+(to) = ^'^{to-i/S), = and $"(to) = 

$'^(to)- Here, J] represents the Lagrangian density for the scalar field $ with 
source J and is defined in terms of the Lagrangian L as 

J] = V^L[<I>] + M . 

It is clear that this generating functional corresponds to a path integral along the 



complex time contour of Fig. 3.1 



The Lagrangian density may now be broken up into pieces quadratic in the field 
and into interaction terms, Ci^t, which are of higher order in the field variables. The 
quadratic integrals may be evaluated explicitly, leaving the following expression in 



terms of the currents and a set of Green's functions[75 



y r ' 

Z[J+, J-, Jf^] = exp I i y" d^x y* c/t{A„,(-i5/(5J+) - £i„,(-i(5/(5J")} 



dtCU-i5/5J^) 

to 

X exp!^^Jj'xid^X2r{xi)G''\xuX2)J\x2)^ . (3.9) 

Here, the subscript c in the last integrals indicates that the time portion of the 
integrals are to be taken over the contour in Fig. p.l| and the superscripts a and b 
have values given by +, — , and f3 corresponding to the location on this contour of 
the times ti and t2- 

This expression simplifies somewhat for the computation of the real time cor- 
relation functions in which we will be interested. For such computations, the con- 
tribution from terms in which a = b = (3 in ( |3.9[ ) cancels between the numerator 
and denominator of Eq. ( p.5| ). In addition, in the limit to ~^ "OO the contributions 
from terms in which either a = (3 oib = (3 vanish for real time correlation functions 
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as a consequence of the Riemann-Lebesgue lemma. We arrive at the generating 
functional 1761 



X exp J^d^xid^X2r{xi)G''\xi,X2)J\x2)^ , (3.10) 

where a and b now take on the values + and — . It should be noted that the finite 
temperature initial condition on p will enter through the boundary conditions on 
the Green's functions as we will see below. 

The Green's functions which appear along the contour may be written in terms 
of the Wightman function G{xi,X2) = ($(xi)$(x2)). They are 

G++(Xi,X2) = G(xi,X2)e(ti-t2)+G'(x2,xi)0(t2-ti) , (3.11) 

G—{Xi,X2) = G{x2,Xi)e{ti-t2) +G{Xi,X2)Q{t2-ti) , (3.12) 

G^^{xi,X2) = G~'^{xi,X2) = —G{xi,X2) ■ (3.13) 

These Green's functions are subject to the periodicity condition in imaginary time: 

G{x2, t2] xi, to) = G(fi, to - il3; X2, ^2) • (3-14) 

The Wightman function G{xi,X2) is constructed from the homogeneous solutions 
of the quadratic form, subject to the periodicity conditions of Eq. ( p.l4|) and 



appropriate initial conditions. We will explicitly construct these solutions within 
the the Hartree and large N approximations for a scalar field below. In addition, 
we will derive the equations of motion for the field and its expectation value. An 
important feature of the Green's functions G"* of Eqs. ( |3.11| )-( p^.l3| ) is that at equal 
times, ti = t2, they are all equivalent and equal to G{xi, X2). 
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3.2 Field Equations in FRW 

We now derive the equations of motion for both a real scalar field and an 0{N) 
invariant scalar field in the vector representation. As was discussed in the in- 
troductory chapter, it is often necessary to study inflationary dynamics using a 
non-perturbative framework. To this end, we derive the equations of motion for 
the real and 0{N) fields within the self-consistent Hartree[^ and the large 
limit respectively. 

These approximations have a variety of important features which make them 
useful in the treatment of inflation. 

1. Within the realm of validity of perturbation theory, the approximations re- 
produce the results of a one loop perturbative expansion. 

2. The approximations remain self-consistent even when amplitudes become 
non-perturbatively large (order 1/A). 

3. As we shall see, the approximations are renormalizable for potentials con- 
taining terms up to fourth order in the scalar field. 

4. The equations of motion preserve the property of covariant conservation of 
the energy-momentum tensor. 

5. The approximations result in local integro-differential evolution equations 
which are quite amenable to integration on a computer. 

In addition to these properties, the large A^ approximation preserves the Ward 
identities corresponding to the 0{N) symmetry. As a consequence, in the case of 
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spontaneous symmetry breaking, the late time dynamics will satisfy Goldstone's 
theorem. 

We will assume, here and throughout, the spatially flat, homogeneous and 
isotropic expanding Robertson- Walker cosmology described in Ch. 2. 

3.2.1 Real Scalar Field: Hartree Dynamics 



We now derive the equations of motion for a real scalar field with Lagrangian (|2.13| ) 



within the self-consistent Hartree approximation!^^. We will find it convenient to 



break the field $ into its expectation value plus a fluctuation about this value: 

$(f,t) = (j){t) +tlj{x,t) , (3.15) 
= (<f(f,t)). (3.16) 

Here, depends only on time due to space translation invariance as is consistent 
with the metric (|2.4|). By definition {i/j{x,t)) = 0; the expectation values in these 
expressions are defined by Eq. ( |3.4| ). 

The Hartree approximation consists of replacing i/j'^'^ by ci(?/'^)"'~^'?/'^ + C2{ip'^)'^ 
and by c?j{ip'^)'^ip^ where the constant factors Cj will be determined below. In 

this way, the Lagrangian ( |2.13| ) is rendered quadratic in the field variable ip, and 
the functional integrals in ( ^.8| ) become purely Gaussian. 

To determine the coefficient Ci, we first make the replacement 

where the numerical factor appearing before the first term on the right hand side is 
determined by matching the n-point functions of if). Next, we use Wick's theorem 
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to write 



from which we find 



/ 2"-i(^-l)! ' 

(2n)! 
~ 2"(n- 1)! ■ 

Using Wick's theorem for the quantity {ip^"-) then gives us the value of C2: 

(2n)!(n - 1) 



C2 



2"n! 



The coefficient C3 is determined in the same way and is found to be given by 

(2n+ 1)! 



C3 



2"n! 



The Hartree factorization may be summarized as follows 



1\n - 1)! 



2"n! 



2"n! 



(3.17) 



Given this factorization, any function -F($) becomes 



where we use the notation 



(3.18) 



(3.19) 



The equations of motion for the scalar field are given by the tadpole condition 



(V') = 1 77]. We find the equation for the mean field 



(3.20) 
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Defining the Fourier transform of the Wightman function by the expression 



G{x, t;x,t) 



(27r)= 



where we have again used the property of space translation invariance. The quan- 
tity (ip'^lt)) = —iG^^{x,x) is constructed from the mode functions obeying the 
equation 

[Gk{t,t)r'fk{t) = 0, (3.21) 



together with the boundary conditions ( |3.14|) , where the operator [Gk] ^ is given 
by the quadratic form appearing in the generating functional (|3.8|) . Explicitly, the 
fk{t) obey 



fk{t) = 



(3.22) 



To this equation we will have to append an appropriately chosen set of initial 
conditions for the mode functions fk{t)- 

Implementing the periodicity condition ( p.l4| ), we find 

i 1 



Gk{t, t 



2uOk 1 - e-^^fe 



fk{t)rk{t) + e-p-'^rk{t)fk{t) 



(3.23) 



and we finally arrive at the expression 

{^■\t)) • 



2(27r) 



.|/,(t)pcoth 



(3.24) 



where T is the initial temperature. Note that in the T limit, the hyperbolic 
cotangent has the value 1. The frequency Uk is given by 

ul = e + a\t,) {V^^ (0(to) + ^(f , to))) - . (3.25) 
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We have included in this definition of the initial frequency a term proportional 
to the Ricci scalar. This term arises when one considers initial conditions corre- 
sponding to the adiabatic vacuum state in conformal time||29|] as was discussed in 



Sec. 2.4. Using such a vacuum state for the mode functions leads to the following 
set of initial conditions on the fk{t)'- 

1 



fkito) 



A(^o) = (-""^^-luj^m^) . (3.26) 

This completes our set of equations of motion of the matter fields within the 
Hartree approximation. We will return to the specific case of the theory in 



Sec. 3.3. 



3.2.2 0{N) Scalar Field: Large N Dynamics 

We now consider $ to be an component scalar field in the vector representation 
of the 0{N) symmetry group. The Lagrangian is 

L[$] = ^V^^x) ■ Vf,${x) - V {${x) ■ $(x)) , (3.27) 

Again, we will separate the mean field value from the fiuctuation. Without loss 
of generality, we choose the expectation value of the i^^ component of the field 
= VNa6i^i such that the mean value of all components of $ except the first 
vanish. We therefore decompose the field as follows: 

$(f, t) = (ViVCT(t) + a{x, t), 7f{x, t)) . (3.28) 

The quantity a is the fiuctuation in the "sigma" direction, whereas vr is the fiuctua- 
tion in the A^— 1 "pion" directions. Since each of the pion directions are equivalent, 

54 



we may write 



7f(x,t) =7r(x,t)(l,l,---,l) . 

For simplicity, we will write the 0{N) generating functional keeping only terms 
on the + time contour in Eq. ( |3.8|) and we will drop the contour superscripts. When 
we explicitly construct the Green's functions below, it will be important to recall 
that these are closed time path Green's functions. 

Defining the currents J and J corresponding to the a and vf fields respectively, 
the generating functional becomes 



Z[J, J] 



-{VNa + cr)n{y/Na + a) - Ja 



+ -jr . njf - Jn + NV i— 



(3.29) 



with □ = V^V^. We now make the functional pion integral Gaussian by intro- 
ducing an auxiliary field x = tt'^/N. The potential term in ( p.29|) may be written 
as 



i / d*xNv(^ 



(3.30) 



It will be convenient to impose the delta function constraint of (|3.3CI|) by introduc- 
ing a Lagrange multiplier field A. The generating functional now reads 



Z[J,J] 



VcrVnVxVX exp 



2Ja + 7f ■ DTf - 2J ■ TT 



X exp <; 2 / d'^xX ( X - ^ 



iN / d^xV 



d X 



[\^a + a)a(\^a + a) 



Na + 

'n 



X 



(3.31) 
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Since we are interested in the dynamics to leading order in the expansion, 
we now Taylor expand the potential V about the quantity a/^/N, keeping terms to 
order The resulting generating functional is quadratic in the vf and a fields. 

Letting /3 = X/N, defining the functions 



G-'lx] ^ n + 2V'{a' + x)+2a'V"{a' + x), 
Fj[x] = J-VNa(n + 2V'{cr^ + x)] 



(3.32) 
(3.33) 
(3.34) 



and performing the Gaussian integrals, we arrive at the functional 

Z[J, J^] = J VxVPexp |z(iV - 1) [^Tr log G-i[/5] + ^J.G[/3] 
X exp |i 

X exp liN I (fx (3x - V{a'^ + x) 



'-Tr\ogG-'[x] + \Fj[x]G,[x]Fj[x] 



(3.35) 

Since each of the pion directions is equivalent, we have set each of the — 1 
components of the current J equal to .]■„. 

Taking to be large, we can now evaluate the x ^i^nd j3 integrals in ( |3.35| ) via 
the saddle point method. Varying the exponentials with respect to {3 and x, we 
find the respective stationary point values x ^-nd j3: 



X 



-iG[P] - .J^G[P]G[P]J^, 



—Fj[x]GMFj[x] 



(3.36) 
(3.37) 



where the 0{1/N) contribution is dropped in the latter variation in determining 
p. 
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For consistency, we drop the 0{1/N) terms to arrive at the generating func- 
tional for normahzed n-point functions 



W[J,U = -lnZ[J,J,] 



+^Fj[x]G,[x]Fj[x] . 



(3.38) 



We now impose the trivial tadpole relation (a) = 0, which yields the equation of 
motion for the zero mode a: 



6J 



^ \d + 2V' (a^ + X) 



j=o 



a . 



(3.39) 



Prom this, we see that we may replace Fj[x] — J and, because the corresponding 
contribution to iy[J, J^] and to ^ is of (l/iV) with respect to the leading order 
contributions, the current J decouples completely from the leading order dynamics. 
To leading order, the generating functional W[J, J^] becomes 

1 



W[J„] = iN 



Tr\ogG-'\p] + -JM^]J^ + / d^^M - J d^xV{a^ + x) 



The two point function is given by 



(tt^) = i-i) 



.,2 6'W[J^] 



-iG[Pi 



(3.40) 



(3.41) 



J^=o 



where, setting = 0, 



P = V'ia' + x), 



(3.42) 
(3.43) 
(3.44) 
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The field a obeys the equation of motion given by Eq. ( p.39| ). 

We now recall that the Green's function, ( |3.32| ), is defined on the OTP contour 
of fig. ( |3.1| ) and that we must impose the periodicity condition ( |3.14| ). The results 
of this procedure [see Eq. ( |3.23D ] are as follows 



U + 2V' (a\t) + {Tx\t)))]a{t)=Q 



2(27r 



^l/.WPcothg 



(3.45) 
(3.46) 



where Uk is defined by 

col = e + 2a\to)V' {^\to) + {n^to)) 
and the mode functions obey the equation 



a2(^o)7^(^o) 



(3.47) 



D + 2V'{a\t) + {7r\t))) Mt)=0. 



(3.48) 



To this set of equations of motion we add the initial conditions on the mode 
functions, which again take the form ( p.26| ) with the frequencies Uk now defined 
by Eq. (13:471) . 

3.2.3 Density Matrix Formalism 

An alternative method of computing the equations of motion without explicit 
computation of the Green's functions is to make an ansatz for the form of the 
density matrix and compute its time evolution via the Liouville equation (|3.1|) . 
Expectation values of operators may then be evaluated explicitly using Eq. (|3.4]) . 
To ensure that such a computation is tractable, it is necessary to choose a Gaussian 
form for the density matrix. We will find that this choice is consistent with the 
approximations we will use to study the dynamics. 
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We restrict our analysis to a real scalar field. The case of an A^-component field 
is completely analogous. With the metric ( p.4| ) and Lagrangian density ( |2.13| ), the 
momentum operator conjugate to is defined as 



n(x) 



(3.49) 



The Hamiltonian may therefore be written as 



li^ix) ait) ,\2 q,,,,,^, 

-A^ + ^(V$(x)) +a\t)V{^x)) 



(3.50) 



_2a?{t) ' 2 

In the Schrodinger picture, the canonical momentum operator at an arbitrary time 
to has the representation 

5 



n(x,to) = -i- 



(3.51) 



5^x,to) ' 

As before, we decompose the field $(x) into its expectation value (f){t) and a 
fluctuation term ipi^x) as in Eq. (|3.16| ). In terms of the density matrix p(t). 



<P{t) = {<^{x,t)) =Tr[pm{x)] 



Again, we have by definition 



(z^(f,t)) = 0. 



(3.52) 



(3.53) 



In terms of these new variables, the Hamiltonian becomes 

1 52 a{t) 



Hit) 



d'^x 



2a3(t) 5^/>2(x) 



W^p{x] 



+ a\t) 



V m)) + V' (0(t)) ^Pix) + -V" (0(t)) ^P\x) + ■ 



(3.54) 



Keeping terms only to quadratic order in iIj{x) corresponds to the first order term 
in a loop expansion. Alternatively, the higher order terms represented by the 
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ellipses of ( |3.54| ) may be approximated as in Eqs. ( |3.17| ), leading to the Hartree 
approximation. We will return to this point shortly. 

It will be convenient to introduce a Fourier representation for the Hamiltonian. 
Letting fl represent the spatial volume, we define 

1 



^{x, t) 



Defining the effective frequencies 



k 



d^xiJi;{t)e-''-^ 



(3.55) 



(3.56) 



the Hamiltonian (|3.54|) becomes, to quadratic order (1-loop), 



H{t) = Qa^t)V{(j){t)) + Y, 



6^ 



+ 



(3.57) 



We now introduce an explicit Gaussian ansatz for the form of the density ma- 



trix: 



P 



Pk{t) 



A;=0 



Uk{t) exp \ -^Afc(t)$fe$_fc - ^A^(t)l>fe$_fc - Bk{t)^k^-k 



(3.58) 



This particular form is dictated by the Hermiticity of the density matrix, 

which also leads to the result that the mixing parameter Bk{t) is real. We choose 
the factorized form since there is no mode-mode mixing in the Gaussian approxi- 
mation. 
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The non-equilibrium evolution of the system is now determined by the Liouville 
equation (p.lD, which reads 



dp ^,i)]t 



dt 



E 



2a3(t) 
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+ a\t)v' m))^h,,{i'k-i'-%) 



+ 



a\t). 



Wl{t) - i,^^i,_A p ^, ^; t . (3.59) 



Inserting our explicit ansatz for the density matrix yields the equations of motion 
for the parameters: 



Mk{t) 
Mk{t) 

ifc(t) 



—I 



^k ~ Bl 
a\t) 

Bk{t) 



a\t)Wl{t) 



Bk{t) = -tJ^{Ak{t)-Al{t))^2 



Bk{t)Aki{t) 



(3.60) 
(3.61) 
(3.62) 



Here, we have introduced the notations AkR and Aki for the real and imaginary 
parts of Ak respectively. 

There are two important notes to make regarding these equations: 

• An initial pure state with Bk{to) = remains pure. 



The time evolution is unitary as indicated by the relations 



d ( Bk{t) 



d 



dt \AkR{t) 



0, 



dt \JAkR{t)+Bk{t)^ 



. 



(3.63) 



(3.64) 
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The relation ( |3.63| ) suggests that we introduce a new function 



Aki = Aki(t) , 



(3.65) 
(3.66) 



where the constants Ck and depend on the choice of initial conditions. These 
definitions yield the following Ricatti-type equation of motion for the quantity 
Ak = AkR + iAki'- 



Mi) 



Am 



This equation may be linearized by a change of variables. We define 

Ak{t) = -ia\t)- 



fkit) 



(3.67) 



(3.68) 



The mode functions fk(t) then obey the equation: 

m+^^^fkit) + wmfkit) = o. 



(3.69) 



To again impose the Hartree approximation, we make a Hartree factorization of 
the effective frequencies by writing 



a?{t) 



+ {v (0(t) + v^(f,t))) 



where the latter term on the right hand side is factorized as in Eq. ( p.l^ ) . In terms 
of these mode functions, the original parameters of the density matrix are given 
by: 

1 



AkR{t) 

Aki{t) 
Bk{t) 



Ck\fk{t)\'' ' 



mm 
1 

DklfkW ■ 



(3.70) 
(3.71) 
(3.72) 
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The two-point fluctuation becomes 



1 



(27r)3 2(Afc^ + 5, 



(27r)3 2 Ck + D, 



(3.73) 



and the equation of motion of the expectation value of $ becomes 



0(t)+3^0(t) + V (0(t)) + 



mt)) = . 



(3.74) 



It is clear from these expressions that for a particular choice of the constants Ck and 
Dk and the initial conditions on the fk, we have reproduced the set of equations 
of motion determined in the last section using a Green's function approach. In 
particular, the choice 



Ck 
Dk 
/fc(to) 

fkito) 



tanh(/5ci;fc) , 
-sinh(/5c^fc) 



'd{to 



a{t(i 



+ iuJk fk{to) , 



along with the appropriate choice of initial normalization 

1/2 



o?{tQ)uk ^ , / (3ujk 
tanh 



TT 



(3.75) 
(3.76) 
(3.77) 

(3.78) 



(3.79) 



yields a density matrix ( p.58|) which is in local thermodynamic equilibrium at the 
initial time at a temperature T = 1/ (3. The two-point fluctuation becomes: 



i2TiY2{AkR + Bk) 



d'k |/fc(t)p^^^^ fPuk 



{2n 



(3.80) 
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3.2.4 Gravitational Dynamics 

The gravitational sector includes the usual Einstein term in addition to a higher 
order curvature term and a cosmological constant term which are necessary to 
renormalize the theory. The action for the gravitational sector is therefore: 

" n{t) 



Sg = J d'^x Cg = j d^xa\t) 



+ -n\t)-K 



(3.81) 



with K being the cosmological constant (we use K rather than the conventional 
A/SttG to distinguish the cosmological constant from the ultraviolet cutoff A we 
introduce to regularize the theory). In principle, we also need to include the terms 
af3^iy as they are also terms of fourth order in derivatives of 
the metric (fourth adiabatic order), but the variations resulting from these terms 
turn out not to be independent of that of 71^ in the flat FRW cosmology we are 
considering. 

The variation of the action S = Sg + Sm with respect to the metric g^^, gives 
us Einstein's equation 

G,. 



+ aH^, + Kg^, = -T^, , (3.82) 



where Guv is the Einstein tensor given by the variation of ^—glZ^ Hn^ is the higher 



order curvature term given by the variation of y/—gTl'^, and Tn^ is the contribution 



from the matter Lagrangian. With the metric ( p.4|) , the various components of the 
curvature tensors in terms of the scale factor are: 

G-o = -3(a/a)^ (3.83) 

g; = -,^ = -6(^ + ^), (3.84) 

64 



H'o = -6(^11 + -^11-^11'], (3.85) 



a 12 
a, 
a 



Ri^ = -Q[n + 3-7^j . (3.86) 
However, it should be clear that even this extended Einstein's equation ( 3.82| ) 



is not self-consistent in the present case, since the energy momentum tensor on the 
right hand side is a quantum quantity, while the gravitational tensors on the left 
hand side are classical. This leads us to one of the major problems of the physics of 
the early universe: the reconciliation of quantum field theory with the cosmological 
theory based on Einstein's equation of general relativity. In particular, there has 
yet to be constructed a full quantum theory of gravity. One must therefore rely 
on approximations. A promising candidate is semi-classical gravity, based on the 
equation 

^^'^ ' aH,, + Kg,, = -{T,,) , (3.87) 



SttGn 

where the expectation value of the quantum energy momentum tensor acts as a 
classical source to the gravitational field. This equation may be considered the first 
approximation to a more complete expansion which includes expectation values of 
terms of higher order in T^^,. The hope is that, since each additional factor of 
Tfj^iy requires a coefficient which by dimensional analysis is of order 1/Mpi, for 
energy densities sufficiently removed from the Planck scale, Eq. (|3.87|) is a good 
approximation to the complete theory. Of course, without a complete quantum 
theory of gravity with which this approximation could be compared, there is no 



definitive way to judge its vahdity (see the discussion in p3|). 

Eventually, when we have fully renormalized the theory, we will simplify the 
theory by setting an to zero, and we will keep as our only contribution to Kn a 
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piece related to the matter fields which we shall incorporate into T^^. Note that 
the dimensional analysis indicates that a is of order l/Mpi so that there may be 
some justification to neglecting this contribution in determining the gravitational 
equations of motion for scales significantly removed from the Planck scale. 

3.3 The A<^^ Theory 

We now specialize the equations of motion we have derived above for very general 
potentials to those appropriate to the theory. 

It is convenient to treat the Hartree and large N approximations simultane- 
ously. We may do so by writing the potential as 

Vmx)) = K+^{m' + mt)) ^'{x) + , (3.88) 

where in the Hartree approximation, the quantity 7 has the value 1/3 and we 
take = 1, while in the large N limit, 7 = 1 and it is to be understood that 
$2 ^ The constant quantity K is a, cosmological constant contribution which 
does not directly affect the equations of motion of $(x) but will be necessary in 
the treatment of spontaneous symmetry breaking. 

Given this potential and using a unified notation by letting — > and (tt^) — > 
{ijj'^) in the equations of motion for the 0{N) scalar field, we arrive at the following 
set of equations of motion in the theory: 

(3.89) 

fk{t) = 0, 

(3.90) 
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mt)) = /^i^wi'-^h(^) . (3.91) 



The initial conditions on the mode functions are 



fk{to) 



1 



fk{to 



with the frequencies uJk{to) given by 



a{to 



- iuJk{to) 



fkito) 



(3.92) 



^^k(to) 



1/2 



(3.93) 



Notice that the only difference between the expressions for the case of a real 
scalar field in the Hartree approximation and the case of an A^-component field 
in the large N limit is a factor of three appearing in the self interaction term in 
the equations for the zero mode. However, as we will see, the two approximations 
are describing theories with distinct symmetries and there will be qualitative dif- 
ferences in the results. An important point to note in the large equations of 



motion is that the form of the equation for the zero mode (|3.89| ) is the same as 
for the A; = mode function ( p.90| ). It will be this identity that allows solutions of 
these equations in a symmetry broken scenario to satisfy Goldstone's theorem. 

3.4 Renormalization and Energy Momentum 

Upon examination of the equal time correlator ( |3.91| ), one finds that the integral 
is divergent and thus must be regulated. This can be done by a number of meth- 
ods, but we require a scheme which is amenable to numerical calculation. We 
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therefore introduce a large momentum cutoff which renders the integral finite, and 
one finds that it is possible to remove the terms depending both quadratically 
and logarithmically on the cutoff by a renormalization of the parameters of the 
theory |]26||. After these terms are subtracted, the cutoff may be taken to infinity, 
and the remaining quantity is both physical and finite. A similar process is re- 
quired to regulate the expressions for the energy density and the pressure as will 
be described in more detail below [p^] . 

In terms of the variables introduced in section II above, the renormalization of 
( p.91|) proceeds almost identically in the Hartree and large approximations. We 
begin by undertaking a WKB analysis of the mode functions fk{t) to reveal their 
large k behavior. We introduce the WKB modes 

, . A „ , . 1 



to / y'^k 

We then have 

fk = {l-^k)Dl + ^kDk. 

The Dk satisfy the same equation of motion as the fk, whereas the function Rk 
obeys the Riccati equation 

Rk{t) + Rlit) + S^Rkit) + ^ + M\t) = . (3.94) 
a[t) a'^[t) 

Writing Rk{t) in the form 

„ / N k „ , , a(t) „ , , a'^(t) „ , , a^(t) „ , , a'^(t) „ , , 

Rkit) = -t^ + Roit) - i^Ri{t) + -^R2{t) - i^R3it) + -^RA{t) + ■ ■ • , 

we may determine Rk{t) term by term by matching coefficients of different powers 
of k in the Riccati equation (|3.94|) . In this way, we find 

m) = (3.95) 
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R2{t) 



1 
2 



6 



i?3W+44j|^3(t) + 2i?i(t)i?(t) 



(3.96) 
(3.97) 
(3.98) 
(3.99) 



The coefficient 7^ is determined from the initial conditions on the mode func- 
tions ( |3.92| ). We find for the real and imaginary parts respectively 



7fc,« = 1 + Oil/k^) , 7fc,, = 0(1/P) . 

Using these results, we find the high k behavior of the mode functions and their 
derivatives. In terms of the effective mass term for the large limit and defining 
the quantities 



B{t) = a\t)[M\t)-n/6] 



■^B j2(.\ , -^3/12/ 



(3.100) 

(3.101) 



where the subscript B stands for bare quantities, we find the following large k 
behavior for the case of an arbitrary scale factor a{t) (with a(0) = 1): 

1 1 



+ 



ka^{t) 2k^a^{t) 



B{t) 



d 



ok^am fm'+<^)jM'^m] + oii/k^) 

5(') + 0(l/P) , 



(3.102) 



k 1 

+ 



B{t) + 2a' 



a^{t) 2ka\t) 
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(3.103) 



fk{t)m+h{t)m 



a{t) 



ka^{t)a{t) Ak^a^it) 



(3.104) 



Because there are no counterterms to utilize in the renormahzation of the the- 
ory, the renormahzation conditions on the mass, couphng to the Ricci scalar and 
coupling constant are obtained from the requirement that the frequencies that 
appear in the mode equations are finitel2^, i.e: 



A 



-m], + iBn{t) + '-^cl>^{t) + '-^{i,'\t))B = -mi + iRn{t)+'-^cf\t) + '-^{ij\t))R 



X 



B /„/,2 



X 



« /„/,2/ 



We find the following renormahzation scheme: 

Xn A' 



TUB + 
Xb = 

is = iR + 



167r2a2(t) ""^ 
Xr 

1 -7Aijln(A/K)/167r2' 
Xb 



1 + ^ ln(A//t) 



IGvr' 



,(^^-1/6) ln(A/ft:), 



(3.105) 

(3.106) 
(3.107) 
(3.108) 



where k is the renormahzation point. The subtracted equal time correlator is now 
given by 



(27r) 



2ka^(t) 



(3.109) 



Notice in particular that in order for the renormahzation of the mass to be time 
independent, we must require that the cutoff A be fixed in physical coordinates and 
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therefore have the form A oc a{t). The abihty to make a renormahzation which is 
independent of time is non-triviaL In particular, choosing initial conditions in with 
respect to the comoving Hamiltonian instead of the initial conditions in conformal 



time (|3.92|) we use here, results in renormalizations which depend explicitly on the 



initial state [^. Further discussions of the significance of the choice of initial state 
may be found in [^. It is worth noting that for a{t) = 1 the renormalizations 
of the mass, quartic coupling, and correlator are the same as those in Minkowski 
space. 

Our treatment of the renormahzation of the energy momentum tensor is similar 



to the approach of [^, extended to the non-perturbative Hartree and large N 
approximations. The expressions for the expectation values of the energy density. 



e, and the trace of the stress energy, e — 3p, where p is the pressure density are[p9[]: 
N 2^ 2 ^ 8 ^ 27A 0^ ^ ''a^^ 

- ^G',{i;') + 6^-{^i;), (3.110) 

a 

'-j^ = + 2mV + + ^ - + (# + + 3^<^ 

- (1 - 60(^') + W)') + (2 - 60m2(^2) - eGj:(l - 60(^') 
+ ^[(2-60^'{ij') + {^Y-6^{^'){i;')n], (3.111) 

where we have used the equations of motion in deriving this expression for the 
trace ( pTTTl) . The quantities and G[] are given by Eqs. and (^^), (ip^) 

is given by Eq. (|3.91| ), {iP'^)r by p.l09| ), and we have defined the following integrals: 

(wf) = /^^v.wr, (3.112) 
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■\m\'. 



(3.113) 



2(27r)3' 

The composite operator (ipip) is symmetrized by removing a normal ordering con- 
stant to yield 



d^k d\fk{t)\' 



(3.114) 



2'^""' ' 4 J (27r)3 dt ■ 

Each of these integrals is divergent and must be regularized. We proceed in 
the same manner as above, imposing an ultraviolet cutoff, A, and computing the 
high k expansions of the fk, this time to fourth order in 1/k. We find the following 
divergences in e and e — 3p: 



A' 



+ 



A2 



div 



+ 



\yi{K/k) 



167r2 



ml 



2{iR-l/Q)Gl + ml + 



X 



^1 j2 



+ 2{in - 1/6)G° [ml + ^{<f + (^2)^)) + {in - 



R) 



(3.115) 



e — 



3p\ 



A2 



div 



IGvr^a^ 



2(e/j-l/6)Gi: + 12(ei?-l/6)- 



+ 2m^ + Ah(0' + (^')ii) 



+ 



ln(A//t) 
167r2 



A 



- 6(eii-l/6) 



■^^(0^ + (^%)+3^^|(0^ + (^^).) 



(3.116) 



The energy momentum is made finite by subtraction of the divergent pieces 
( |3.115 ) and ( 3.116| ) from the expressions for the energy density ( ^.11(J| ) and the 
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trace ( |3.111| ). Within the context of covariant regularization schemes such as 



dimensional regularization and covariant point splitting|^, such a procedure has 
been shown to adequately renormalize couplings appearing in the semi-classical 
Einstein's equation, Eq. ( |3.87D . 



The renormalizations of Newton's constant, the higher order curvature cou- 
pling, and the cosmological constant are given by the condition of finiteness of the 
semi-classical Einstein equation: 

^0 ^^^H',+KB9'o + {T!^)B = ^ + cynH',+Kn9'o + {T!^)R- (3.117) 



Finally we arrive at the following set of additional renormalizations pO 

1 _ 1 / 1\ 1\ 2 ln(A/K) 

" S^^iVG^ " H^"" " 6 j 16^ + l^^ " 6 j ""^ ' 

(3.118) 



As expected, the logarithmic terms are consistent with the renormalizations found 



using dimensional regularization |^ Again, we set aji = and choose the 



renormalized cosmological constant such that the vacuum energy is zero in the 
true vacuum. We emphasize that while the regulator we have chosen does not 
respect the covariance of the theory, the renormalized energy momentum tensor 
defined in this way nevertheless retains the property of covariant conservation in 
the limit when the cutoff is taken to infinity. 

In the numerical analysis of inflation, the logarithmic subtractions can be ne- 
glected because of the coupling A < 10^^^. Using the Planck scale as the cutoff 
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and the inflaton mass ttir as a renormalization point, these terms are of order 
Xlii[Mpi/mR] < 10~^°, for m > 10^ GeV . An equivalent statement is that for 
these values of the coupling and inflaton masses, the Landau pole is well beyond 
the physical cutoff Mpi. Our relative error in the numerical analysis is of order 
10~^, so that our numerical study is insensitive to the logarithmic corrections. 
Though these corrections are fundamentally important, they are small enough to 
be negligible in our numerical analyses. 
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Chapter 4 



New Inflation 



4.1 Introduction and Motivation 



With the introduction of the first model of inflation by Guth|T^ came an immediate 
problem. This 'old' inflation model failed in that it did not include a mechanism 
to gracefully exit from the inflationary phase. The basic problem was that 'old' 
inflation relied on a first order phase transition in which the inflaton begins in a 
meta-stable inflationary state and then tunnels to move to the true minimum of 
the potential resulting in a radiation or matter dominated phase, but the rate for 
such a phase transition to take place was too small. 

Soon, it was realized that inflation could occur in a second order phase transi- 
tion as well, and a model based on the Coleman- Weinberg theory was built which 
had no exit problem since the filed evolved continuously to the true vacuum. This 
'new' inflation model[^, however, had a problem of its own. Since the coupling 



in the Coleman- Weinberg theory is a gauge coupling of order one, the metric per- 
turbations produced in this model were much too large to be consistent with the 
observed homogeneity of the CBR[^. Nevertheless, the basic idea of inflation 
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Figure 4.1: A typical new inflation potential. The field $ begins near the top of 
the central hill. 

proceeding via a second order transition is sound and a number of models, such as 
the model with negative mass squared, are successful in producing inflation 
with all the desired properties. 

In this chapter, we examine the non-equilibrium dynamics of a new infla- 
tionary quantum field theory. The expectation value is initially near the symmetric 
point of the potential (see Fig. the disordered state, and the fluctuations about 
this value are assumed to be small relative to the vacuum expectation value (vev), 
V ~ n/VX (where n = \m\). 

Because this initial state is in a region of the potential with negative curva- 
ture, the field will grow in time due to the presence of spinodal instabilities ||41|| . 
This growth continues until the fluctuations become of the same order as the vev, 
i.e. ($^(x, t)) ~ t>^ ~ /^^/A. Until this time, when the field has become non- 
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perturbatively large, inflation occurs. 

On the initial state: dynamics of phase transitions 

The issue of initial conditions on the inflaton is of particular interest in treatments 



of new inflation. In much of the literature P, ^ , the inflaton is treated as a classical 
field (f)ci{t) which is taken to be initially 'near' the origin. Further, the uncertainty 
relation of quantum mechanics is sometimes invoked to yield a requirement that 

Mto) > H/2TT. 

How this classical field relates to the quantum field t), however, has been 
somewhat unclear. The straightforward identification of the classical field (j)d{t) 
with the expectation value = t)) is clearly not correct, as there are 

no restrictions on the initial value 0(to) and, in particular, a natural value in a 
phase transition 0(to) = would yield the incorrect result of an unending period 
of inflation if (f){t) were identifled with a classical inflaton fleld. 

One is led to the following questions: 

• Assuming it can be done, how does one construct the classical inflaton from 
the underlying quantum fleld theory? 

• What, if anything, does this construction tell one about the initial state of 
the effective classical theory? 

As we shall see, the answers to these questions lie in the quantum fluctua- 
tions ip{x,t) of the fleld. Let us take the example of a symmetric initial condition 
(pito) = The zero mode remains zero for all time (pit) = 0. However, examina- 
tion of Eq. (|3.90| ) for negative mass squared indicates a spinodal instability for low 
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k which results in growth of the non-zero momentum modes which make up the 
quantum fluctuations {ip'^it)), Eq. ( |3.91| ). Therefore the dynamics becomes fluc- 
tuation driven with (ip'^lt)) eventually becoming non-perturbatively large (that is, 
(ip^lt)) ~ /i^/A). One may then hope that the classical picture of new inflation may 
be constructed from the quantum evolution by the identification of the classical 
field with some component of the quantum fluctuations. 
However, questions remain: 

• How does the gravitational background respond to non-perturbatively large 
quantum fluctuations? 

• How can one extract (small?) metric perturbations from non-perturbatively 
large field fluctuations? 

Remaining still are the details of how the classical field might be constructed from 
the quantum field. In this chapter, we address each of these questions within a 
fully consistent quantum field theory context. 

Before we move on to our analysis of new inflation, there is one more point 
regarding initial conditions that deserves attention. The standard picture of a 
second order phase transition in the early universe has the field $(x, t) in thermal 
equilibrium at some temperature above critical, T > T^, such that the initial state 
is one of positive mass squared. This localizes the field about the origin. Then, as 
the finite temperature contributions to the effective mass are redshifted away, the 
phase transition takes place and the field may evolve toward its asymptotic state. 

As we shall see, however, the self-coupling of the inflaton is constrained to be 
A < 10"^^ by the magnitude of the temperature fluctuations of the GBR. It has 
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therefore been argued that the field is too weakly interacting to be in thermal 
equilibrium above the critical temperature. It is therefore not clear whether there 
is any 'natural' initial state for the field in new infiation. 

As discussed in Sec. 3.3, we will simply assume an initial density matrix describ- 
ing either a thermal state at temperature T = or a vacuum state of the initial 
Hamiltonian (an initial state in thermal equilibrium at a temperature T = 0). The 
assumption of an initial equilibrium vacuum state are essentially the same used 
by Linde|^0[, Vilenkin[pl[], as well as by Guth and Pi[|l^ in their analyses of the 



quantum mechanics of inflation in a fixed de Sitter background. 

A final note is that the initial frequencies ujk ( |3.92| ) are imaginary for low 



k modes in the vacuum state and need to be modified. This may be done in 
a variety of ways with little effect on results[^, |3^. Here, we choose a smooth 
interpolation between low-k modes with modified frequencies and the high k modes 
which remain in the conformal vacuum state with frequencies ujk 



.^ + A^^(to)tanh^ |A^^(to)| 



(4.1) 



4.2 Analysis 

We will work in the large N approximation throughout this chapter. Except for 
the very late time dynamics with 0(to) 7^ 0, the evolution in the Hartree case is 
almost identical. For 0(to) = 0, there is no difference between the Hartree and 
large N dynamics. 

It is convenient to introduce the following dimensionless quantities and defini- 
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tions, with /i = \mfi\, 

, H k Wk , 

r = ^^t ; h = - ; g=- u, = — ; , = — , (4.2) 

^'(r) = ^0^(t) ; 9m = ^mt))n ; /,(r) = . (4.3) 

Choosing = (minimal couphng) and the renormahzation point k = |mij| 
and setting a(0) = 1, the equations of motion become: 



— 3A— 



l + r/2 + (^S(. 



r/(r) 



0, 
0, 



/,(0) = [-/i(0)-z^,]/,(0), 



where 



Sfr) 



q^dq 



a r 



6m| _ 



(4.4) 
(4.5) 

(4.6) 



2q-^ \ m)j 

The initial conditions for 77 (r) will be specified later. An important point to 
notice is that the equation of motion for the g = mode coincides with that of 
the zero mode ([4.4| ). Furthermore, for ?7(r 00) 7^ 0, a stationary (equilibrium) 
solution of the Eq. 



is obtained when the sum rule 25, 27, 29 



- 1 + r/2(cx)) + ^S(cx)) = 



(4.7) 



is fulfilled. This sum rule is nothing but a consequence of Goldstone's theorem and 
is a result of the fact that the large approximation satisfies the Ward identities 
associated with the 0{N) symmetry, since the term — 1 + r/^ + gT, is seen to be 
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the effective mass of the modes transverse to the symmetry breaking direction, i.e. 
the Goldstone modes in the broken symmetry phase. 

In terms of the zero mode rjlr) and the quantum mode function given by 
Eq. ( [4. 51) we find that the Friedmann equation for the dynamics of the scale factor 
in dimensionless variables is given by 



, 9 / X ,9 / N ,9 4:7lNm% 

h\r) = Ahl enir) ; /^^ = (4.^ 
and the renormahzed energy and pressure are given by: 



enir) = + \ {-1 + ii' + gi:)' + 



2 ' 4 
2 

,2 



Jq'dq [|//-5«(g,r) 
+ ^(|//-e(g-l)5(^)(g,r); 
ip + eUr) = ^-^[f,' + glq'dq[\f,\'-S^'\q,r) 



(4.9) 



+ £^{\f,\'-QiQ-l)S^'\q,r)) |, (4.10) 

where the subtractions S^^^ and S^"^^ are given by the right hand sides of Eqs. ( p.l03| ) 
and ( p.l02| ) respectively. 

The renormalized energy and pressure are covariantly conserved: 

en{T) + 3h{T){p + e)R{T) = 0. (4.11) 

In order to provide the full solution we now must provide the values of r]{0), 
fj{0), and Hq. Assuming that the inflationary epoch is associated with a phase 
transition at the GUT scale, this requires that A^m^/A/j ~ (10^^ GeV )^ and 
assuming the bound on the scalar self-coupling A/j ~ 10^^^ — 10^^^ (this will be 



seen later to be a compatible requirement), we find that ho ^ A^^/^ which we will 
take to be reasonably given by /iq ~ 1 — 10 (for example in popular GUT's ^ 20 
depending on particular representations). 

We will begin by studying the case of most interest from the point of view of 
describing the phase transition: rj^O) = and ?7(0) = 0, which are the initial con- 
ditions that led to puzzling questions. With these initial conditions, the evolution 
equation for the zero mode Eq. ( |4.4| ) determines that //(r) = by symmetry. 

4.2.1 Early time dynamics 

Before engaging in the numerical study, it proves illuminating to obtain an esti- 
mate of the relevant time scales and an intuitive idea of the main features of the 
dynamics. Because the coupling is so weak {g ~ 10^^^ <ti 1) and because after 
renormalization the contribution from the quantum fluctuations to the equations 
of motion is finite, we can neglect all the terms proportional to g in Eqs. ( [4.9|) and 

(PD- 

The evolution equations for the mode functions are those for an inverted oscil- 



lator in de Sitter space-time, which have been studied by Guth and Pi[|T^]. One 
obtains the approximate solution 



h{t) ^ ho, 



/9 1 
4 + /I!' 



(4.12) 



where J±y{z) are Bessel functions, and Aq and Bq are determined by the initial 
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conditions on the mode functions: 
1 nq 



/uJ^2hQ sinz/TT 



tUJ, 



g — ^ho ^ f q 



ho 



J' I r- 



IT q 



/un2hQ sin z/tt 



iujq - I hp J /q^ 
q \ho 



ho 



(4.13) 
(4.14) 



After the physical wavevectors cross the horizon, i.e. when qe '^"'^//i-o ^ 1 we 
find that the mode functions factorize: 



fgir) 



B„ 



2hn 



(4.15) 



Til-u) \ q J 

This result reveals a very important feature: because of the negative mass squared 
term in the matter Lagrangian leading to symmetry breaking (and u > 3/2), we 
see that all of the mode functions grow exponentially after horizon crossing (for 
positive mass squared u < 3/2, and they would decrease exponentially after horizon 
crossing). This exponential growth is a consequence of the spinodal instabilities 
which is a hallmark of the process of phase separation that occurs to complete 
the phase transition. We note, in addition that the time dependence is exactly 
given by that of the q = mode, i.e. the zero mode, which is a consequence 
of the redshifting of the wavevectors and the fact that after horizon crossing the 
contribution of the term q^/a^{T) in the equations of motion become negligible. 
We clearly see that the quantum fluctuations grow exponentially and they will 
begin to be of the order of the tree level terms in the equations of motion when 
5fE(r) 1. At large times 

S(r)^^2^/lo)/^^e(2^-^)'^«^ 

with J-'{ho) a finite constant that depends on the initial conditions and is found 
numerically to be of 0{1) [see Fig. 4.11|. 
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In terms of the initial dimensionful variables, the condition 5'S(r) ^ 1 translates 
to {ip"^ {x , t)) ^ 2m|./Aij, i.e. the quantum fluctuations sample the minima of the 
(renormalized) tree level potential. We find that the time at which the contribution 
of the quantum fluctuations becomes of the same order as the tree level terms is 



estimated to be||27 
1 



(2z/-3)/io 



In 



ghlT^iho) 



-ho In 



ghlJ^^{ho 



+ O{l/ho). (4.16) 



At this time, the contribution of the quantum fluctuations makes the back reaction 
very important and, as will be seen numerically, this translates into the fact that 
Tg also determines the end of the De Sitter era and the end of inflation. The total 
number of e-folds during the stage of exponential expansion of the scale factor 
(constant ho) is given by 



2z/ 



In 



9 hi T^ih) 



-hi In 
2 ° 



9 hi J'^h) 



+ 0(1) 



(4.17) 



For large ho we see that the number of e-folds scales as hi as well as with the 
logarithm of the inverse coupling. These results (|4.15| ) - ( [4.171 ) will be confirmed 
numerically below and will be of paramount importance for the interpretation of 
the main consequences of the dynamical evolution. 

4.2.2 Classical or quantum behavior? 

Above we have analyzed the situation when 77(0) = (or in dimensionful variables 
0(0) = 0). The typical analysis of inflaton dynamics in the literature involves 
the classical evolution of (f){t) with an initial condition in which 0(0) is very close 
to zero (i.e. the top of the potential hill) in the 'slow-roll' regime, for which 
^ 3H(f). Thus, it is important to quantify the initial conditions on 0(t) for 
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which the dynamics will be determined by the classical evolution of and those 
for which the quantum fluctuations dominate the dynamics. We can provide a 
criterion to separate classical from quantum dynamics by analyzing the relevant 
time scales, estimated by neglecting non-linearities and backreaction effects. We 
consider the evolution of the zero mode in terms of dimensionless variables, and 
choose r]{0) 7^ and f]{0) = 0. (?7(0) 7^ simply corresponds to a shift in origin of 
time). We assume ^^(O)^ << 1 which is the relevant case where spinodal instabilities 
are important. We find 

7]{t) ^ 7]{0) e(^-i)^«" . (4.18) 

The non-linearities will become important and eventually terminate inflation 
when ?7(r) ^ 1. This corresponds to a time scale given by 

r. « . (4.10) 

[i^ - f ) ho 

If Tc is much smaller than the spinodal time given by Eq. ([4.16|) then the classical 
evolution of the zero mode will dominate the dynamics and the quantum fluctua- 
tions will not become very large, although they will still undergo spinodal growth. 
On the other hand, if Tc ^ Tg the quantum fluctuations will grow to be very large 
well before the zero mode reaches the non-linear regime. In this case the dynamics 
will be determined completely by the quantum fluctuations. Then the criterion for 
the classical or quantum dynamics is given by 

77(0) ^ y/g =^ classical dynamics 

77(0) ^ ^/g ho =^ quantum dynamics (4.20) 

or in terms of dimensionful variables 0(0) ^ Hq leads to classical dynamics and 
0(0) <^ Hq leads to quantum dynamics. 

85 



However, even when the classical evolution of the zero mode dominates the 
dynamics, the quantum fluctuations grow exponentially after horizon crossing un- 
less the value of (f){t) is very close to the minimum of the tree level potential. In 
the large N approximation the spinodal line, that is the values of (j){t) for which 
there are spinodal instabilities, reaches all the way to the minimum of the tree 
level potential as can be seen from the equations of motion for the mode functions. 
Therefore even in the classical case one must understand how to deal with quantum 
fluctuations that grow after horizon crossing. 

4.2.3 Numerics 

The time evolution is carried out by means of a fourth order Runge-Kutta routine 
with adaptive step sizing while the momentum integrals are carried out using an 11- 
point Newton-Cotes integrator. The relative errors in both the differential equation 
and the integration are of order 10~®. We flnd that the energy is covariantly 
conserved throughout the evolution to better than a part in a thousand. Figures 



44- [4.4| show 5'S(r) vs. r, /i(r) vs. r and ln(|/g(r)| ) vs. r for several values of 
q with larger q's corresponding to successively lower curves. Figures [4.5| and [4.6| 
show p{T)/e{T) and the horizon size h~^{T) for g = 10^^^ ; ?7(0) = ; ?7(0) = 
and we have chosen the representative value ho = 2.0. 

Figures [4.2| and |4.3| show clearly that when the contribution of the quantum 
fluctuations 5'S(r) becomes of order 1 inflation ends, and the time scale for 5'S(r) 
to reach 0{1) is very well described by the estimate ([4.16|) . From Fig. |4.2| we see 



that this happens for r = ~ 90, leading to a number of e-folds A^e 180 which 
is correctly estimated by ( [4.1(j| , [4.17| ). 
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Figure 4.4: \n{\fg{T)\'^) vs. r, for 77(0) = 0,77(0) ^ 0,g ^ lO'^^ho = 2.0 for 
q — 0.0, 5, 10, 15, 20 with smaller q corresponding to larger values of ln{\ fq{T)\^). 
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Figure 4.5: p/e vs. r, for 77(0) = 0, 77(0) = 0, = 10-l^ ho = 2.0. 
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Figure [4.4| shows clearly the factorization of the modes after they cross the 



horizon as described by Eq. ( |4.15| ). The slopes of all the curves after they become 
straight lines in Fig. ff^ is given exactly by (2z/ — 3), whereas the intercept depends 



on the initial condition on the mode function and the larger the value of q the 
smaller the intercept because the amplitude of the mode function is smaller initially. 
Although the intercept depends on the initial conditions on the long-wavelength 
modes, the slope is independent of the value of q and is the same as what would 
be obtained in the linear approximation for the square of the zero mode at times 
long enough that the decaying solution can be neglected but short enough that 
the effect of the non-linearities is very small. Notice from the figure that when 
inflation ends and the non-linearities become important all of the modes effectively 
saturate. This is also what one would expect from the solution of the zero mode: 
exponential growth in early-intermediate times (neglecting the decaying solution), 
with a growth exponent given by (z/ — 3/2) and an asymptotic behavior of small 
oscillations around the equilibrium position, which for the zero mode is r/ = 1, but 
for the 7^ modes depends on the initial conditions. All of the mode functions 
have this behavior once they cross the horizon. We have also studied the phases 
of the mode functions and we found that they freeze after horizon crossing in the 
sense that they become independent of time. This is natural since both the real 
and imaginary parts of fqij) obey the same equation but with different boundary 
conditions. After the physical wavelength crosses the horizon, the dynamics is 
insensitive to the value of q for real and imaginary parts and the phases become 
independent of time. Again, this is a consequence of the factorization of the modes. 
The growth of the quantum fluctuations is sufficient to end inflation at a time 
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Figure 4.6: l//i(r) vs. r, for r/(0) = 0,?7(0) = 0,^ = 10'^^, ho = 2.0. 



given by Tg in Eq. (|4.16| ). Furthermore Fig. |4.5| shows that during the inflationary 
epoch p{t)/s{t) ^ —1 and the end of inflation is rather sharp at r<j with p{T)/e{T) 
oscillating between ±1 with zero average over the cycles, resulting in matter dom- 
ination. Figure [4.6| shows this feature very clearly; /i(r) is constant during the 
De Sitter epoch and becomes matter dominated after the end of inflation with 
h-\T) ^ |(r - r,). There are small oscillations around this value because both 
p{t) and e(r) oscillate. These oscillations are a result of small oscillations of the 
mode functions after they saturate, and are also a feature of the solution for a zero 
mode. 

All of these features hold for a variety of initial conditions. As an example, we 
show in Figs. [4.7| - [4.1U| the plots corresponding to Figs. |4.2| - f4.5| for the case of an 
initial Hubble parameter of ho = 10. 
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Figure 4.7: gJl vs. r, for r/(0) = 0, ?)(0) = 0, c/ = 10"^^ /iq = 10.0. 
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Figure 4.8: H{t) vs. t, for 77(0) = 0, 77(0) = 0, ^ = 10"^^, ho = 10.0 
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Figure 4.9: ln(|/g(r)|2) vs. r, for r/(0) = 0,?7(0) = 0,g = W-^^,ho = 10.0 for 
q = 0.0, 5, 10, 15, 20 with smaller q corresponding to larger values of ln(|/g(r)p). 
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Figure 4.10: p/e vs. r, for r/(0) = 0,?7(0) = 0, c/ = 10-^^ /iq = 10.0. 



92 



4.3 Zero Mode Assembly 

This remarkable feature of factorization of the mode functions after horizon cross- 
ing can be elegantly summarized as 

fk{t)\k^h{t)^Ho = 9{q, ho)Mr), (4.21) 

with kph{t) — k e~^°* being the physical momentum, g{q, ho) a complex constant, 
and /o(t) a real function of time that satisfies the mode equation with q = and 
real initial conditions which will be inferred later. Since the factor g{q, ho) depends 
solely on the initial conditions on the mode functions, it turns out that for two 
mode functions corresponding to momenta ki, k2 that have crossed the horizon at 
times ^1 > t2, the ratio of the two mode functions at time t, {ts > t > ti > t-^) is 



Then if we consider the contribution of these modes to the renormalized quan- 
tum fluctuations a long time after the beginning of inflation (so as to neglect the 
decaying solutions), we find that 

where 'small' stands for the contribution of mode functions associated with mo- 
menta that have not yet crossed the horizon at time r, which give a perturbatively 
small (of order A) contribution. We find that several e-folds after the beginning 
of inflation but well before inflation ends, this factorization of superhorizon modes 
implies the following: 
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gjq'dq\fl{r)\ ^ |Co|Vo'(r), (4.22) 
gjq^dq\fl{r)\ ^ \C,\'fl{r), (4.23) 
^ I - ^^o^r), (4.24) 

where we have neglected the weak time dependence arising from the perturbatively 
small contributions of the short-wavelength modes that have not yet crossed the 
horizon, and the integrals above are to be understood as the fully renormalized 
(subtracted), finite integrals. For = 0, we note that ( [4.22| ) and the fact that 
/o('r) obeys the equation of motion for the mode with g = leads at once to the 
conclusion that in this regime [(^^(r)]^ = |Co|/o(t) obeys the zero mode equation 
of motion 

^ + 3h^-l + i\Co\foir)Y |Co|/o(r) = . (4.25) 

It is clear that several e-folds after the beginning of inflation, we can define an 
effective zero mode as 

Veffi^) = 9^i^)^ oi" ill dimensionful variables, (peffit) = {'^'^{^■,t))R ^ (4.26) 

Although this identification seems natural, we emphasize that it is by no means 
a trivial or ad-hoc statement. There are several important features that allow an 
unambiguous identification: i) [{ip'^{x^ '^))r\ is a fully renormalized operator product 
and hence finite, ii) because of the factorization of the superhorizon modes that 
enter in the evaluation of [('?/'^(x, t))/?], 0e//(^) ([4.26|) obeys the equation of motion 
for the zero mode, iii) this identification is valid several e-folds after the beginning 
of inflation, after the transient decaying solutions have died away and the integral 
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in {ilj'^{x,t)) is dominated by the modes with wavevector k that have crossed the 
horizon at t{k) <^ t. Numerically we see that this identification holds throughout 
the dynamics except for a very few e-folds at the beginning of infiation. This 
factorization determines at once the initial conditions of the effective zero mode 
that can be extracted numerically: after the first few e-folds and long before the 
end of infiation we find 

0e//(t) = <f>eff{0) e(^-i)^o* , (4.27) 

where we parameterized 



to make contact with the literature. As is shown in Fig. ^4.11| , we find numerically 
that J-'{Ho/m) ^ 0{1) for a large range of 0.1 < Hq/itl < 50 and that this quantity 
depends on the initial conditions of the long wavelength modes. 
Therefore, in summary, the effective composite zero mode obeys 



^ + 34-l + ,=„(r) 



Veff(T) = ; iW/(t = 0) = (i^ - 2) ,fc„(0) , (4.28) 



where ?7e//(0) = ^^^^ (0) is obtained numerically for a given by fitting the 
intermediate time behavior of 5'S(r) with the growing zero mode solution. 

Furthermore, this analysis shows that in the case = 0, the renormalized energy 
and pressure in this regime in which the renormalized integrals are dominated by 
the superhorizon modes are given by 

, , 1Nm% fin 1 / 9 \2l , , 

- ^{2^'// + 4(-l + ^'//) I' ^4-29) 

{V^e)n - ^ (4.30) 

95 



2.0 • 



1 .5 



1 -O 



0.5 




O.O 



1 O.O 



20.0 



30.0 



40.0 



50.0 



Figure 4.11: J^{H/m) vs. H, where J^{H/m) is defined by the relation (f)eff{0) 
{H/2TT)J^{H/m) [see Eqs. ( ^M ) and (^) l. 



where we have neglected the contribution proportional to l/a^(r) because it is 
effectively red-shifted away after just a few e-folds. We found numerically that 
this term is negligible after the interval of time necessary for the superhorizon 
modes to dominate the contribution to the integrals. Then the dynamics of the 
scale factor is given by 



^eff 



(4.31) 



We have numerically evolved the set of effective equations ( [4.28| , [4.31|) by ex- 
tracting the initial condition for the effective zero mode from the intermediate time 
behavior of gIl{T), and found a remarkable agreement between the evolution of r/g^j 
and 5'S(r) and between the dynamics of the scale factor in terms of the evolution 
of rjefflr), and the full dynamics of the scale factor and quantum fluctuations 
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Figure 4.12: r^^^/r) vs. r, for r/e//(0) = 3.94 x 10-^77e//(0) = 0.317r/e//(0), = 
10~^^, ho = 2.0. The initial conditions were obtained by fitting the intermediate 
time regime of 5'S(r) in fig. ([4.21). rjeffir) is the solution of Eq. ( [4.2^ ) with these 
initial conditions. 



within our numerical accuracy. Figs. 4.12 and 4.13 show the evolution of Tjlffir) 
and /i(r) respectively from the classical evolution equations ([4.28|) and ( |4.31|) us- 
ing the initial condition rjeffiO) extracted from the exponential fit of 5'S(r) in the 
intermediate regime. These figures should be compared to Figs. [4.2| and [4.3| . We 
have also numerically compared p/e given solely by the dynamics of the effective 
zero mode and it is again numerically indistinguishable from that obtained with 
the full evolution of the mode functions. 

This is one of the main results of our work. In summary: the modes that 
become superhorizon sized and grow through the spinodal instabilities assemble 
themselves into an effective composite zero mode a few e-folds after the beginning 
of inflation. This effective zero mode drives the dynamics of the FRW scale factor, 
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Figure 4.13: /i(r) vs. r, obtained from the solution of eqns. ( [4.28| ) and ( |4.31| ) with 
the conditions of fig. (^4.12|) ). 



terminating inflation when the non-linearities become important. In terms of the 
underlying fluctuations, the spinodal growth of superhorizon modes gives a non- 
perturbatively large contribution to the energy momentum tensor that drives the 
dynamics of the scale factor. Inflation terminates when the root mean square 
fluctuation probes the equilibrium minima of the tree level potential. 

This phenomenon of zero mode assembly, i.e. the 'classicalization' of quantum 
mechanical fluctuations that grow after horizon crossing is very similar to the in- 



terpretation of 'decoherence without decoherence' of Starobinsky and Polarski||82 



The extension of this analysis to the case for which ri{0) 7^ is straightforward. 



Since both ri{T) and y gT,{T) = |Co|/o(t) obey the equation for the zero mode, 
Eq. ( [4.4| ), it is clear that we can generalize our definition of the effective zero mode 
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to be 



Veff{r) = ^r/2(r) +^S(r) . 
which obeys the equation of motion of a classical zero mode: 

Veffir) = 



dr 



(4.32) 



(4.33) 



If this effective zero mode is to drive the FRW expansion, then the additional 
condition 



• 2 



?77o -2w/o/o + r/7o =0 



(4.34) 



must also be satisfied, which is simply the condition that the contribution of the 
kinetic term in the effective theory is the same as that in the full theory. One can 
easily show that this relation is indeed satisfied if the mode functions factorize as in 
( [4.21| ) and if the integrals ( ^.22| ) - ( |4.24| ) are dominated by the contributions of the 
superhorizon mode functions. This leads to the conclusion that the gravitational 
dynamics is given by eqns. ( f4.29| ) - ( [4. 311 ) with ?7e//(T) defined by ( |4.32| ). 

We see that in all cases, the full large quantum dynamics in these models 
of inflationary phase transitions is well approximated by the equivalent dynam- 
ics of a homogeneous, classical scalar field with initial conditions on the effective 
field ?7e//(0) > ^/ghoJ-'{hQ). We have verified these results numerically for the 
field and scale factor dynamics, finding that the effective classical dynamics repro- 
duces the results of the full dynamics to within our numerical accuracy. We have 
also checked numerically that the estimate for the classical to quantum crossover 
given by Eq. ( [4.2(]| ) is quantitatively correct. Thus in the classical case in which 
r]{0) ^ VX ho we find that rjeffij) = rj^r), whereas in the opposite, quantum case 
flefAr) -- 



r . 
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This remarkable feature of zero mode assembly of long-wavelength, spinodally 
unstable modes is a consequence of the presence of the horizon. It also explains 
why, despite the fact that asymptotically the fluctuations sample the broken sym- 
metry state, the equation of state is that of matter. Since the excitations in the 
broken symmetry state are massless Goldstone bosons one would expect radiation 
domination. However, the assembly phenomenon, i.e. the redshifting of the wave 
vectors, makes these modes behave exactly like zero momentum modes that give 
an equation of state of matter (upon averaging over the small oscillations around 
the minimum). 

Sub-horizon modes at the end of inflation with q > 6^°'^" do not participate 
in the zero mode assembly. The behavior of such modes do depend on q after the 
end of inflation. Notice that these modes have extremely large comoving q since 
li^^hoTs y 10^^. As discussed in Ref. p9[ such modes decrease with time after 



inflation as ~ l/a(r). 

Making sense of the fluctuations 

Having recognized the effective classical variable that can be interpreted as the 
component of the field that drives the FRW background and rolls down the classi- 
cal potential hill, we want to recognize unambiguously the small fluctuations. We 
have argued above that after horizon crossing, all of the mode functions evolve 
proportionally to the zero mode, and the question arises: which modes are assem- 
bled into the effective zero mode whose dynamics drives the evolution of the FRW 
scale factor and which modes are treated as perturbations? In principle every 
k mode provides some spatial inhomogeneity, and assembling these into an 
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effective homogeneous zero mode seems in principle to do away with the very inho- 
mogeneities that one wants to study. However, scales of cosmological importance 
today first crossed the horizon during the last 60 or so e-folds of inflation. Recently 
Grishchuk |]83|l has argued that the sensitivity of the measurements of AT/T probe 



inhomogeneities on scales ~ 500 times the size of the present horizon. Therefore 
scales that are larger than these and that have first crossed the horizon much ear- 
lier than the last 60 e-folds of inflation are unobservable today and can be treated 
as an effective homogeneous component, whereas the scales that can be probed 
experimentally via the CMB inhomogeneities today must be treated separately as 
part of the inhomogeneous perturbations of the CMB. 

Thus a consistent description of the dynamics in terms of an effective zero mode 
plus 'small' quantum fluctuations can be given provided the following requirements 
are met: a) the total number of e-folds Ne ^ 60, b) all the modes that have crossed 
the horizon before the last 60-65 e-folds are assembled into an effective classical 
zero mode via (peffit) = [<Po{t) + {'ip'^{x, t))^]^ , c) the modes that cross the horizon 
during the last 60-65 e-folds are accounted as 'small' perturbations. The reason 
for the requirement a) is that in the separation (f){x,t) = (f)gff{t) + 6(f){x,t) one 
requires that 6(j){x,t)/(j)eff{t) <ti 1. As argued above, after the modes cross the 
horizon, the ratio of amplitudes of the mode functions remains constant and given 
by e^'^-V^^ with AN being the number of e-folds between the crossing of the 
smaller k and the crossing of the larger k. Then for 6(j){x, t) to be much smaller than 
the effective zero mode, it must be that the Fourier components of 5(f) correspond 
to very large /c's at the beginning of inflation, so that the effective zero mode 
can grow for a long time before the components of 50 begin to grow under the 
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spinodal instabilities. In fact requirement a) is not very severe; in Figs. [4.2| -^l6 



we have taken ho = 2.0 which is a very moderate value and yet for A = 10"^^ 
the inflationary stage lasts for well over 100 e-folds, and as argued above, the 
larger for fixed A, the longer is the inflationary stage. Therefore under this 
set of conditions, the classical dynamics of the effective zero mode (f)eff{t) drives 
the FRW background, whereas the inhomogeneous fluctuations 6(f){x,t), which are 
made up of Fourier components with wavelengths that are much smaller than the 
horizon at the beginning of inflation and that cross the horizon during the last 60 
e-folds, provide the inhomogeneities that seed density perturbations. 

4.4 Scalar Metric Perturbations 

Having identifled the effective zero mode and the 'small perturbations', we are now 
in position to provide an estimate for the amplitude and spectrum of scalar metric 
perturbations. We use the formulation by Mukhanov, Feldman and Branden- 
berger[^, a short review of which was provided in Sec. 2.5. 

We are interested in determining the dynamics of $fc for those wavevectors that 
cross the horizon during the last 60 e-folds before the end of inflation. During the 
inflationary stage the numerical analysis yields, to a very good approximation 

Hit) ^ Ho ; 0e//(t) = 0e//(O) e^'^-i)^^*, (4.35) 

where Ho is the value of the Hubble constant during inflation, leading to 



(u-2)Hot 



Ho J ^ \ Ho ^ 
The coefficients and bk are determined by the initial conditions. 
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P = u-1. (4.36) 



Since we are interested in the wavevectors that cross the horizon during the 
last 60 e-folds, the consistency for the zero mode assembly and the interpretation 
of 'small perturbations' requires that there must be many e-folds before the last 
60. We are then considering wavevectors that were deep inside the horizon at the 
onset of inflation. Mukhanov et. al.|^ show that is related to the canonical 

'velocity field' that determines scalar perturbations of the metric and which is 
quantized with Bunch-Davies initial conditions for the large fc-mode functions. 
The relation between and v and the initial conditions on v lead at once to a 



determination of the coefficients and bk for k >> Hq 



ak 



Sir 



(0), 



TT 1 



2Hr,k 







(4.37) 



Thus we find that the amplitude of scalar metric perturbations after horizon 
crossing is given by 



Mt)\=k^Mt)\^\ 



(0) 



2Ho 
k 



(4.38) 



The time dependence of \Sk(t)\ displays the unstable growth associated with the 
spinodal instabilities of superhorizon modes and is a hallmark of the phase transi- 



tion. This time dependence can be understood from the constraint equation (2.42) 



that relates the Bardeen potential to the gauge invariant field fiuctuations|P5 
which in longitudinal gauge are identified with 6(f){x,t). 

Since the right hand side of ( p.42| ) is proportional to (pejf/Mpi ^ 1 during the 
inflationary epoch in this model, we can neglect the terms proportional to and 



on the left hand side of (|2.43|) , in which case the equation for the gauge invariant 



scalar field fluctuation is the same as for the mode functions. In fact, since 



A5« i 



IS 
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gauge invariant we can evaluate it in the longitudinal gauge wherein it is identified 
with the mode functions fk(t). Then, absorbing a constant of integration in the 
initial conditions for the Bardeen variable, we find 

" i^''^ ^"'^'"> *' t^-^") 



and using that oc e('^~3/'^)^ot g^^^ j^^^g^f, after horizon crossing fk{t) oc e^'^"^/^)^"*, 
one obtains at once the time dependence of the Bardeen variable after horizon 
crossing. In particular the time dependence is found to be oc e('^'^-^)Hot ^ jg ti^en 
clear that the time dependence is a refiection of the spinodal (unstable) growth of 
the superhorizon field fiuctuations. 

To obtain the amphtude and spectrum of density perturbations at second hori- 
zon crossing we use the conservation law, Eq. ( 2.44| ), which is valid after horizon 



crossing of the mode with wavevector k. Although this conservation law is an exact 
statement of superhorizon mode solutions of Eq. ( |2.41| ), we have obtained solutions 



assuming that during the infiationary stage H is constant and have neglected the 
H term in Eq. ( |2.41| ). Since during the infiationary stage, 

= -W, " "° (^^)' "° 

and (f)/(j) ^ Hq, the above approximation is justified. We then see that (plfjit) oc 
g(2i/-3)//ot -^viiiQj^ ig same time dependence as that of $fe(t). Thus the term 
proportional to 1/(1 +p/e) in Eq. (|2.44]) is indeed constant in time after horizon 



crossing. On the other hand, the term that does not have this denominator evolves 
in time but is of order (l+p/e) = —2H/3H^ ^ 1 with respect to the constant term 
and therefore can be neglected. Thus, we confirm that the variable ^ is conserved 
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up to the small term proportional to (1 +p/e)(^k which is negligible during the 
inflationary stage. This small time dependence is consistent with the fact that we 
neglected the H term in the equation of motion for 

However, upon second horizon crossing it is straightforward to see that ^k{tf) ~ 
0. The reason for this assertion can be seen as follows: Eq. ( p.43| ) shows that at 



long times, when the effective zero mode is oscillating around the minimum of 
the potential with a very small amplitude and when the time dependence of the 
fluctuations has saturated (see Fig. $fc will redshift as ~ ^i^nd its 

derivative becomes extremely small. 

Using this conservation law, assuming matter domination at second horizon 
crossing, and $fc(i/) ~ 0|^, we find 



12r(z.)v/^ f2HoY-- 

where T{HQ/m) determines the initial amplitude of the effective zero mode ( |4.27| ). 
We can now read the power spectrum per logarithmic k interval 

Vs{k) = 141' oc A;-2('^-t) = fc'^^-i . (4.42) 

leading to the index for scalar density perturbations 

= 1 - 2 (^z/ - ^) . (4.43) 

For HQ/m ^ 1, we can expand u — 3/2 as a series in m? / Hq in Eq. (|4.41|) . 
Given that the comoving wavenumber of the mode which crosses the horizon n 
e-folds before the end of inflation is = ifoe*-^"""-' where Ne is given by (|4.17|) , 
we arrive at the following expression for the amplitude of fluctuations on the scale 
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corresponding to n in terms of the De Sitter Hubble constant and the couphng A: 



{AAA) 

Here, 7 is Euler's constant. Note the exphcit dependence of the amphtude of 
density perturbations on V^- For n ~ 60, the factor exp(nm^/3ifo) is (9(100) 
for Hq/itl = 2, while it is 0{1) for Hq/iti > A. Notice that for Hq/tti large, the 
amplitude increases approximately as {HQ/mY, which will place strong restrictions 
on A in such models. 

We remark that we have not included the small corrections to the dynamics of 
the effective zero mode and the scale factor arising from the non-linearities. We 
have found numerically that these nonlinearities are only significant for the modes 
that cross about 60 e-folds before the end of inflation for values of the Hubble 
parameter Ho/rriR > 5. The effect of these non-linearities in the large limit 
is to slow somewhat the exponential growth of these modes, with the result of 
shifting the power spectrum closer to an exact Harrison-Zeldovich spectrum with 
Hg = 1. Since for Ho/mji > 5 the power spectrum given by ( [4.43|) differs from 



one by at most a few percent, the effects of the non-linearities are expected to 
be observationally unimportant. The spectrum given by ( [4.41|) is similar to that 



obtained in references pl|, although the amplitude differs from that obtained 
there. In addition, we do not assume slow roll for which (z/ — |) <^ 1, although 
this would be the case if Ne ^ 60. 

We emphasize an important feature of the spectrum: it has more power at 
long wavelengths because u — 3/2 > 0. This is recognized to be a consequence of 
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the spinodal instabilities that result in the growth of long wavelength modes and 
therefore in more power for these modes. This seems to be a robust prediction of 
new inflationary scenarios in which the potential has negative second derivative in 
the region of field space that produces inflation. 

It is at this stage that we recognize the consistency of our approach for separat- 
ing the composite effective zero mode from the small fluctuations. We have argued 
above that many more than 60 e-folds are required for consistency, and that the 
small fluctuations correspond to those modes that cross the horizon during the last 
60 e-folds of the inflationary stage. For these modes Ho/k = e~^"**(^) where t*{k) 
is the time since the beginning of inflation of horizon crossing of the mode with 
wavevector k. The scale that corresponds to the Hubble radius today Aq = 27[/kQ 
is the first to cross during the last 60 or so e-folds before the end of inflation. 
Smaller scales today will correspond io k > k^ at the onset of inflation since they 
will cross the first horizon later and therefore will reenter earlier. The bound on 
\5ka\ oc AT/T < 10^^ on these scales provides a lower bound on the number of 
e-folds required for these type of models to be consistent: 



where we have written the total number of e-folds as Ne = Hq t*{ko) + 60. This in 
turn can be translated into a bound on the coupling constant using the estimate 
given by Eq. (gTT^. 

Since the COBE satellite observes microwave background temperature variation 
with a resolution of a single degree or larger, which is approximately the size of 
the particle horizon at the time of recombination, we may directly compare the 




(4.45) 
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power spectrum for thermal fluctuations measured by COBE with the spectrum 
computed here for fluctuations at second horizon crossing [see Eq. ( p.40| )]. The 
four year COBE DMR Sky Map[^ gives n 1.2 ± 0.3 thus providing an upper 



bound on v 



< z/ - ^ < 0.05 (4.46) 



corresponding to > 2.6. We then find that these values of ho and A ^ 
10~^^ — 10~^^ provide sufficient e-folds to satisfy the constraint for scalar density 
perturbations. 

4.5 Decoherence: Quantum to Classical Transi- 
tion During Inflation 

An important aspect of cosmological perturbations is that they are of quantum 
origin but eventually they become classical as they are responsible for the small 
classical metric perturbations. This quantum to classical crossover is associated 



with a decoherence process and has received much attention p2|, q4 . 

Recent work on decoherence focused on the description of the evolution of 
the density matrix for a free scalar massless field that represents the "velocity 
field" [4^] associated with scalar density perturbations 1 32 1. In this section we study 



the quantum to classical transition of superhorizon modes for the Bardeen variable 
by relating these to the field mode functions and analyzing the full time evolution 
of the density matrix of the matter field. This is accomplished with the identifi- 



cation given by equation ( |4.39| ) which relates the mode functions of the Bardeen 
variable with those of the scalar field. This relation establishes that in the models 
under consideration the classicality of the Bardeen variable is determined by the 
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classicality of the scalar field modes. 

In the situation under consideration, long-wavelength field modes become spin- 
odally unstable and grow exponentially after horizon crossing. The factorization 
( [4. 151) results in the phases of these modes "freezing out". This feature and the 
growth in amplitude entail that these modes become classical. The relation (|4.39| ) 



in turn implies that these features also apply to the superhorizon modes of the 
Bardeen potential. 

Therefore we can address the quantum to classical transition of the Bardeen 
variable (gravitational potential) by analyzing the evolution of the density matrix 
for the matter field. 

To make contact with previous workfS^, Q we choose to study the evolution of 
the field density matrix in conformal time (see Sec. 2.4), although the same features 
will be displayed in comoving time. Note that we use r] to denote conformal time 
throughout this section; it should not be confused with the dimensionless, comoving 
time zero mode fjit). 

The conformal time Hamiltonian operator, which is the generator of transla- 
tions in ?7, is given by 

H, = /rf^x{ln|+i(Vx)^ + V(x)}, (4.47) 

with 11^ being the canonical momentum conjugate to X; = x'. Separating the 
zero mode of the field x 

X{x, v) = Xo{v)ki + X{x, ??), (4.48) 

and performing the large factorization on the fluctuations we find that the 
Hamiltonian becomes linear plus quadratic in the fluctuations, and similar to a 
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Minkowski space-time Hamiltonian with a rj dependent mass term given by 



m' 



(4.49) 



We can now follow the steps and use the results of ref.|2^ for the conformal 
time evolution of the density matrix by setting a{t) = 1 in the proper equations of 
that reference and replacing the frequencies by 



(4.50) 



The expectation value in Eq. ( |4.49| ) and that of the energy momentum tensor are 
obtained in this rj evolved density matrix. [As is clear, we obtain in this way the 
self-consistent dynamics in the curved cosmological background ( p^.29 )]. 

The time evolution of the kernels in the density matrix (see |^) is determined 
by the mode functions that obey 



The Wronskian of these mode functions 



(4.51) 



W{F, F* 



FlFt - F,F': 



(4.52) 



is a constant. It is natural to impose initial conditions such that at the initial r] 
the density matrix describes a pure state which is the instantaneous ground state 
of the Hamiltonian at this initial time. This implies that the initial conditions of 
the mode functions Fkirj) be chosen to be (see 



1 



(4.53) 
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With such initial conditions, the Wronskian ( [4.52| ) takes the value 



W(F, F*) = -2i . 



(4.54) 



The Heisenberg field operators x{x,ri) and their canonical momenta n^(x, 77) 
can now be expanded as: 



(27r)3/2 
d?k 
(27r)3/2 



Ak-x 



ik-x 



(4.55) 
(4.56) 



with the time independent creation and annihilation operators and at obeying 
canonical commutation relations. Since the fiuctuation fields in comoving and con- 
formal time are related by a conformal rescaling given by Eq. ( |2.30| ) it is straight- 
forward to see that the mode functions in comoving time t are related to those in 
conformal time simply as 



fk{t) 



a{ri) ■ 



(4.57) 



Therefore the initial conditions given in Eq. (|4.53|) on the conformal time mode 
functions and the choice a(to) = ^(^o) = 1 imply the initial conditions for the 
mode functions in comoving time given by Eq. (|3.92|) . 
In conformal time, the density matrix takes the form 

k 

- Bkiv) VkC^) ■ V-kiv) + i Hiv) ■ (v-kiv) - V-kiv)) } ; , (4-58) 
Vkiv) = Xkiv) - xoiv) ki ^{k) ; Vkiv) = Xkiv) - xoiv) ki ^{k) ■ 

7fj^(r/) is the Fourier transform of Il^{ri,x). 
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In conformal time quantization and in the Schrodinger representation in which 
the field x is diagonal the conformal time evolution of the density matrix is via the 
conformal time Hamiltonian ( |4.47] ). The evolution equations for the covariances are 
obtained from those given in ref. |^ by setting a{t) = 1 and using the frequencies 
^kiv) = k"^ + Ai'^i'i]). In particular, by setting the covariance of the diagonal 
elements (given by equation (2.20) in see also equation (2.44) of [^]), 



Mv) 



More explicitly 

N'kiv) = J^kim) exp 



AikW) H 



Mk{m] 



\fjkijh) \Fk{rii)\ 



Aikiv) 



A\og\Fk{ri)\ = -ait) -ait) ^log|/,(t)| 



(4.59) 



ARk{r]) 



1 



Wkiv)? a{t)^\Ut)\^ 



(4.60) 



Bk{ri) 



where A^k and Ajk are respectively the real and imaginary parts of Ak and we 
have used the value of the Wronskian ( [4.54| ) in evaluating ( |4.60| ). 

The coefficients Ak{ri) and N'kiv) in the Gaussian density matrix ( [4.5^ ) are 
completely determined by the conformal mode functions Fk{ri) [or alternatively 
the comoving time mode functions fk{t)]- 

Let us study the time behavior of these coefficients. During inflation, a(t) ^ 
e^°^, and the mode functions factorize after horizon crossing. Further, superhorizon 
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modes grow in cosmic time as in Eq. ( [4.151 ): 



where the coefficient can be read from Eq. ( [4.15| ) 



We emphasize that this is a result of the full evolution as displayed from the nu- 



merical solution in Fig. [4.4| . These mode functions encode all of the self-consistent 
and non-perturbative features of the dynamics. This should be contrasted with 
other studies in which typically free field modes in a background metric are used. 



Inserting this expression in Eqs. (|460| ) yields 



Since z/ — | > 1, we see that the imaginary part of the covariance Aikirf) grows 
very fast. Hence, the off-diagonal elements of p[$, $, r^] oscillate wildly after a 
few e-folds of inflation. In particular their contribution to expectation values of 
operators will be washed out. That is, we quickly reach a classical regime where 
only the diagonal part of the density matrix is relevant: 

k 

The real part of the covariance ARk{ri) (as well as any non-zero mixing kernel 
-6fe(?7)[^) decreases as e"^^'^"^-*'^"*. Therefore, characteristic field configurations 
r/g are very large (of order e^'^^^^^^'^). Therefore configurations with field ampli- 
tudes up to 0{e'^'^~'^^^''°^) will have a substantial probability of occurring and being 
represented in the density matrix. 
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Notice that x ~ e*-'^~^)'^o* corresponds to field configurations $ with amphtudes 
of order e('^~i)'*o* [see Eq. ( p.30| )]. It is the fact that z/— | > which in this situation 
is responsible for the "classicalization" , which is seen to be a consequence of the 
spinodal growth of long-wavelength fluctuations. 

The equal-time field correlator is given by 



2(27r)= 

= <i? jl^.m^e''-^'-''^ ■ (4.62) 

and is seen to be dominated by the superhorizon mode functions and to grow as 
^{2v-i)hot^ whereas the field commutators remain fixed showing the emergence of a 
classical behavior. As a result we obtain 

{x{^,v) xi^'yV)) oc a^{t) (peffit) (peffit') G{\x- x'\) + small (4.63) 
where G{\x — x'\) falls off exponentially for distances larger than the horizon P? 



and "small" refers to terms that are smaller in magnitude. This factorization of 
the correlation functions is another indication of classicality. 

Therefore, it is possible to describe the physics by using classical field the- 
ory. More precisely, one can use a classical statistical (or stochastic) field theory 
described by the functional probability distribution ( [4.61| ). 



These results generalize the decoherence treatment given in ref.[£2| for a free 
massless field in pure quantum states to the case of interacting fields with broken 
symmetry. Note that the formal decoherence or classicalization in the density 
matrix appears after the modes with wave vector k become superhorizon sized i.e. 
when the factorization of the mode functions becomes effective. 
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4.6 Conclusions 



It can be argued that the inflationary paradigm as currently understood is one 
of the greatest applications of quantum field theory. The imprint of quantum 
mechanics is everywhere, from the dynamics of the inflaton, to the generation of 
metric perturbations, through to the reheating of the universe. It is clear then that 
we need to understand the quantum mechanics of inflation in as deep a manner as 
possible so as to be able to understand what we are actually testing via the CBR 
temperature anisotropies. 

What we have found in our work is that the quantum mechanics of inflation 
is extremely subtle. We now understand that it involves both non-equilibrium as 
well as non-perturbative dynamics and that what you start from may not be what 
you wind up with at the end. 

In particular, we see now that the correct interpretation of the non-perturbative 
growth of quantum fluctuations via spinodal decomposition is that the background 
zero mode must be redefined through the process of zero mode reassembly. When 
this is done we can interpret inflation in terms of the usual slow-roll approach with 
the now small quantum fluctuations around the redeflned zero mode driving the 
generation of metric perturbations. 

We have studied the non-equilibrium dynamics of a 'new inflation' scenario in 
a self-consistent, non-perturbative framework based on a large N expansion, in- 
cluding the dynamics of the scale factor and backreaction of quantum fluctuations. 
Quantum fluctuations associated with superhorizon modes grow exponentially as a 
result of the spinodal instabilities and contribute to the energy momentum tensor 
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in such a way as to end inflation consistently. 

Analytical and numerical estimates have been provided that establish the re- 
gime of vahdity of the classical approach. We find that these superhorizon modes 
re-assemble into an effective zero mode and unambiguously identify the composite 
field that can be used as an effective expectation value of the inflaton field whose 
classical dynamics drives the evolution of the scale factor. This identification also 
provides the initial condition for this effective zero mode. 

A consistent criterion is provided to extract "small" fluctuations that will con- 
tribute to cosmological perturbations from "large" non-perturbative spinodal fluc- 
tuations. This is an important ingredient for a consistent calculation and interpre- 
tation of cosmological perturbations. This criterion requires that the model must 
provide many more than 60 e-folds to identify the 'small perturbations' that give 
rise to scalar metric (curvature) perturbations. We then use this criterion com- 
bined with the gauge invariant approach to obtain the dynamics of the Bardeen 
variable and the spectrum for scalar perturbations. 

We find that during the inflationary epoch, superhorizon modes of the Bardeen 
potential grow exponentially in time reflecting the spinodal instabilities. These 
long- wavelength instabilities are manifest in the spectrum of scalar density pertur- 
bations and result in an index that is less than one, i.e. a "red" power spectrum, 
providing more power at long wavelength. We argue that this 'red' spectrum is 
a robust feature of potentials that lead to spinodal instabilities in the region in 
field space associated with inflation and can be interpreted as an "imprint" of the 
phase transition on the cosmological background. Tensor perturbations on the 
other hand, are not modified by these features, they have much smaller amplitude 
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and a Harrison- Zeldovich spectrum. 

We made contact with the reconstruction program and vahdated the results 
for these type of models based on the slow-roll assumption, despite the fact that 
our study does not involve such an approximation and is non-perturbative. 

Finally we have studied the quantum to classical crossover and decoherence of 
quantum fluctuations by studying the full evolution of the density matrix, thus 



making contact with the concept of "decoherence without decoherence" p2| which 
is generalized to the interacting case. In the case under consideration decoherence 
and classicalization is a consequence of spinodal growth of superhorizon modes and 
the presence of a horizon. The phases of the mode functions "freeze out" and the 
amplitudes of the superhorizon modes grow exponentially during the inflationary 
stage, again as a result of long-wavelength instabilities. As a result field configura- 
tions with large amplitudes have non-vanishing probabilities to be represented in 
the dynamical density matrix. In the situation considered, the quantum to classical 
crossover of cosmological perturbations is directly related to the "classicalization" 
of superhorizon matter field modes that grow exponentially upon horizon crossing 
during inflation. The diagonal elements of the density matrix in the Schroedinger 
representation can be interpreted as a classical distribution function, whereas the 
off-diagonal elements are strongly suppressed during inflation. 
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Chapter 5 

Explosive Particle Production 
After Inflation 

5.1 Introduction 

The most popular of current models of inflation is probably the chaotic inflation 
models [^. Unlike the old or new inflation models which rely on the inflaton being 
near some local extremum of the potential, typically with 0(to) ^ \^\/ chaotic 
inflation requires initial values of the field which are large, typically with 0(to) ^ 
\m\/\/~\ (see Fig. For such large field values, the potential contribution 

dominates the energy momentum, leading to inflation. 

Unlike the new inflationary scenario examined in the preceding chapter, there 
are no spinodal instabilities in chaotic inflation and the dynamics of the inflation- 
ary phase are relatively straightforward. In particular, during inflation, quantum 
fluctuations remain small and one may simply treat the inflaton as a classical field 

^One may question the validity of such an initial state for the field. The argument in favor of 
chaotic inflation is that one requires only that there exists a single region of the universe of size 
H^^ satisfying the condition at some time t ^ tpi. Since the only (known) restriction on the 
field at the Planck scale is that the energy density satisfies (j) < M^i, such field values exist and 
chaotic inflation occurs at least in some part of the universe. Subsequently, regions which have 
inflated come to dominate the physical volume of the universe. 
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Figure 5.1: A typical chaotic inflation potential. The field $ begins at very high 
field values, far from the minimum. 



identified as the mean value of the quantum field. That is, it is safe to write 
4>d{t) = 4>{t) and be sure that corrections to the classical evolution remain pertur- 
batively small. 

However, as we shall see, the post inflationary evolution may not be so straight- 
forward. 

In addition to other aspects of inflation, the question of the reheating of the uni- 
verse, the transfer of energy from the inflaton to other particle modes, has received 



much attention [Q. The reason for this is that exponential expansion during in- 
flation causes the universe to become very cold, as the energy in flelds other than 
the inflaton is redshifted away. In order to achieve the standard results of nucle- 
osynthesis, and possibly other early processes, the universe must be reheated to 
temperatures above those at which these important processes take place. Early 
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efforts to account for tlie necessary refieating introduced ad hoc decay widths to 
tlie inflaton, assuming tlie energy transfer occurred tlirougli single particle decay 
E7^. However, more recently, it has been realized that there are much more ef- 



ficient processes, those of either spinodal decomposition in the case of a (new) 
inflationary phase transition and parametric amplification in the case of chaotic 



inflationary scenarios p3| |56|, |53, |5g, |60|, |61|, |3g]. Such mechanisms, in 



which the rate of energy transfer grows exponentially, are referred to generally as 
preheating. 

In particular, we will study the process of preheating using a wide range of 
initial conditions. We analyze the dynamics of preheating in fixed radiation and 
matter dominated Friedmann- Robertson- Walker cosmologies both analytically and 
numerically^. 

Our main results are as follows. In the situations we analyze, we find that 
the expansion of the universe allows for significant particle production, although 
this production is somewhat sensitive to the exact expansion rate and is effectively 
shut off for high enough rates. In the case of a symmetry broken potential, we 
determine that in the large N limit, the quantum fluctuations decrease for late 
times as l/a^(t), while, as in Ch. 4, these fluctuations and the zero mode satisfy a 
sum rule consistent with Goldstone's theorem [see Eq. ( |4.7| )]. 

We compute the renormalized energy density, e, and the pressure, p, as a 
function of time. Averaging over the field oscillations, we find immediately after 



^We could easily allow the scale factor to evolve dynamically in the semi-classical approxi- 
mation to gravity, as described in Ch. 3. However, we find this self-consistent evolution to be 
almost identical to the evolution in the appropriate fixed gravitational background. We therefore 
simplify our analysis (and reduce our computational run time) by restricting ourselves to the 
fixed background. 
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preheating a cold matter equation of state {p = 0) in the slow roll scenarios. In 
chaotic scenarios, the equation of state just after preheating is between that of 
radiation [p = e/3) and matter where the matter dominated regime is reached 
only for late times. The time scale over which the equation of state becomes 
matter dominated depends on the distribution of created particles in momentum 
space in addition to the approximation scheme implemented. 

5.2 Evolution 



We focus our study of the evolution on radiation or matter dominated cosmologies, 
as the case for de Sitter expansion has been studied previously[|2^. We write 
the scale factor as a{t) = (t/to)" with n = 1/2 and n = 2/3 corresponding to 
radiation and matter dominated backgrounds respectively. Note that the value of 
to determines the initial Hubble constant since 



a{to) 



n 



We now solve the system of Eqs. ( |3.89| )-( |3.91| ) in the Hartree approximation and 
in the large N limit. We begin by presenting an early time analysis of the slow roll 
scenario. We then undertake a thorough numerical investigation of various cases 
of interest. For the symmetry broken case, we also provide an investigation of the 
late time behavior of the zero mode and the quantum fluctuations. 

It will also be convenient to define the number operator N{t). We do so using 
conformal time quantization, as the particles so defined will most closely match 
the definition of particles in flat spacetime[|29|. We have, at zero temperature. 



Nk{t) 



4 



fkjt) 

fk{to) 



1 + 



h{t) + H{t)h{t) 
fk{t) 



(5.i: 
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Examination of the large k behavior of fk and fk shows that A^^^ is 0{l/k'^) for 
large k, so that the total number of particles 

k 

(SF^^ (5.2) 

is finite. 

In what follows, we will assume minimal coupling to the curvature, = 0. In 
most of the cases of interest, TZ/\m\ <^ 1, so that finite C,r has little effect. 

5.2.1 Early Time Solutions for Slow Roll 

For early times in a slow roll scenario [m^ = —1, r]{To) ^ 1], we can neglect in 
Eq. (p.89|) and in Eq. (|3.90|) both the quadratic and cubic terms in 0(t) as well as 



the quantum fluctuations {'ip'^it))^ [recall that {ip^ito))^ = 0]. Thus, the differential 
equations for the zero mode ( p.89| ) and the mode functions ( p. 901 ) become linear 



equations. In terms of the scaled variables introduced at the beginning of Sec. 4.2, 
with a(r) = (r/ro)" {n = 2/3 for a matter dominated cosmology while n = 1/2 
for a radiation dominated cosmology) we have: 

vir) + -viT)-v{r) = 0, (5.3) 

T 



cP 3n d 
+ ^ + 



dr^ T dr (r/ro)2'^ 
The solutions to the zero mode equation ( |5.3|) are 



/,(r) = 0. (5.4) 



ri{T) = c T-'h{T) + d t-^K,{t) , (5.5) 

where v = {3n — l)/2, and /;^(t) and K^{t) are modified Bessel functions. The 
coefficients, c and d, are determined by the initial conditions on rj. For //(tq) = r/o 
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and ?7(to) = 0, we have: 



^0 ^0 



V 



To 



d 



-m tq 



iu{ro) h{To) 

To 



(5.6) 
(5.7) 



Taking the asymptotic forms of the modified Bessel functions, we find that for 
intermediate times 7]{t) grows as 



9n'^ — 6n 
87 



We see that //(r) grows very quickly in time, and the approximations ( p.3| ) and 
( [5. 41 ) will quickly break down. For the case shown in Fig. (with n = 2/3, 
^(''"o) =0.1, 77(to) = 10"'', and ?7(ro) = 0), we find that this approximation is valid 
up to r — To ~ 10. 

The equations for the mode functions ( ^.4|) can be solved in closed form for the 
modes in the case of a radiation dominated cosmology with n = 1/2. The solutions 
are 



f,{r)=c,e--u(^--'^X^r]+d,e~-M(^--'^,l2r]. 



4 



2 '2' 



2 ' 2' 



(5.9) 



Here, U{-) and M(-) are confluent hypergeometric functions |8^ (in another com- 
mon notation, M(-) = iFi(-)), and the Cq and dg are coefficients determined by 
the initial conditions (|3.92|) on the modes. The solutions can also be written in 
terms of parabolic cylinder functions. 

For large r we have the asymptotic form 



r>l 



dg e"(2r) 



~(3/4+g2ro/2)_ 



TT 



2 F ( 



3 _ qfjo 

4 2 



1 + ( - 



+ Cg e^^ (2r 



,(-3/4+(72r„/2) 



1 + Of- 



(5.10) 
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Again, these expressions only apply for intermediate times before the nonlinearities 
have grown significantly. 

5.2.2 Numerical Analysis 

We now present the numerical analysis of the dynamical evolution of scalar fields in 
time dependent, matter and radiation dominated cosmological backgrounds. We 
use initial values of the Hubble constant such that h{TQ) > 0.1. For expansion rates 
much less than this value the evolution will look similar to Minkowski space, which 



has been studied in great detail elsewhere |^5|, ^ ^ . As will be seen, the equation 



of state found numerically is, in the majority of cases, that of cold matter. We 
therefore use matter dominated expansion for the evolution in much of the analysis 
that follows. While it presents some inconsistency at late times, the evolution in 
radiation dominated universes remains largely unchanged, although there is greater 
initial growth of quantum fluctuations due to the scale factor growing more slowly 
in time. Using the large N and Hartree approximations to study theories with 
continuous and discrete symmetries respectively, we treat three important cases. 
They are 

1. < 0, r/(ro) < 1, 

2. < 0, ■r]{To) > 1, and 

3. > 0, t]{to) > 1. 

In presenting the figures, we have shifted the origin of time such that r ^ r' = 
T — tq. This places the initial time, tq, at the origin. In these shifted coordinates, 
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r + r ^ 



the scale factor is given by 

a(r) = 

where, once again, n = 2/?) and n = 1/2 in matter and radiation dominated 
backgrounds respectively, and the value of f is determined by the Hubble constant 
at the initial time: 

h{To = 0) = ^. 

r 

Case 1: < 0, i]{tq) <^ 1. This is the case of an early universe phase 
transition in which there is little or no biasing in the initial configuration (by biasing 
we mean that the initial conditions break the 1] —r] symmetry). The transition 
occurs from an initial temperature above the critical temperature, T > Tc, which is 
quenched at to to the temperature Tf <^ Tc. This change in temperature due to the 
rapid expansion of the universe is modeled here by an instantaneous change in the 
mass from an initial value = T^/T^^ — 1 to a final value rnj = — 1. We will use 
the value mf = 1 in what follows. This quench approximation is necessary since 
the low momentum frequencies ( p.93|) appearing in our initial conditions ( p.92| ) are 



complex for negative mass squared and small ti{tq). An alternative choice is to use 
initial frequencies ( [4.1D which were used in Ch. 4. We find that different choices 
of initial conditions change the quantitative value of the particle number by a few 
percent, but leave the qualitative results unchanged. 

We plot the the zero mode 77(t), the equal time correlator 5'S(r), the total 
number of produced particles gN{T) [see Eq. ( |5.2| )], the number of particles gNq{T) 
as a function of wavenumber for both intermediate and late times [see Eq. 
and the ratio of the pressure and energy densities p{t)/6{t) (giving the equation 
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of state). 

Figure |5.2| a-e shows these quantities in the large approximation for a matter 
dominated cosmology with an initial condition on the zero mode given by ri{TQ = 
0) = 10"'', ?7(ro = 0) = and for an initial expansion rate of /i(to) = 0.1. This 
choice for the initial value of t] stems from the fact that the quantum fluctuations 
only have time to grow significantly for initial values satisfying tjItq) ^ ^/g; for 
values 77(to) ^ the evolution is essentially classical. This result is clear from 
the intermediate time dependence of the zero mode and the low momentum mode 
functions given by the expressions (|5.8| ) and (|5.10|) respectively. 

After the initial growth of the fluctuation gT, (Fig. ^.2[ b) we see that the zero 



mode (Fig. |5.2| a) approaches the value given by the minimum of the tree level 
potential, f] = 1, while gT, decays for late times as 

C _ C 

^ ~ a2(r) r4/3 " 

For these late times, the Ward identity corresponding to the 0{N) symmetry of 
the field theory is satisfied, enforcing the condition [see (^3)] 

- 1 + r]\T) + gJ:{T) = 0. (5.11) 

Hence, the zero mode approaches the classical minimum as 

C 



7] (t) ~ 1 



Figure [5.2| c depicts the number of particles produced. After an initial burst of 
particle production, the number of particles settles down to a relatively constant 
value. Notice that the number of particles produced is approximately of order 1/g. 
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Figure a 



Figure b 




Figure 5.2: Symmetry broken, slow roll, large N, matter dominated evolution of 
(a) the zero mode r]{T) vs. r, (b) the quantum fluctuation operator 5'S(r) vs. r, 
(c) the number of particles gN{r) vs. r, (d) the particle distribution gNq{T) vs. 
q at T = 149.1 (dashed line) and r = 398.2 (solid line), and (e) the ratio of the 
pressure and energy density p(T)/e(r) vs. r for the parameter values = — 1, 
r/(ro) = 10-^ ?)(ro) = 0, ^ = h{To) = 0.1. 
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In Fig. |5.2| d, we show the number of particles as a function of the wavenumber, 
q. For intermediate times we see the simple structure depicted by the dashed line 
in the figure, while for late times this quantity becomes concentrated more at low 
values of the momentum q. 



Finally, Fig. shows that the field begins with a de Sitter equation of state 
p = —e but evolves quickly to a state dominated by ordinary matter, with an 
equation of state (averaged over the oscillation time scale) p = 0. This last result 
is a bit surprising as one expects from the condition (|5.11| ) that the particles 
produced in the final state are massless Goldstone bosons (pions) which should have 
the equation of state of radiation. However, as shown in Fig. ^^2|d, the produced 
particles are of low momentum, g -C 1, and while the effective mass of the particles 
is zero to very high accuracy when averaged over the oscillation time scale, the 
effective mass makes small oscillations about zero so that the dispersion relation 
for these particles differs from that of radiation. In addition, since the produced 
particles have little energy, the contribution to the energy density from the zero 
mode, which contributes to a cold matter equation of state, remains significant. 



In Fig. |5.3| a-e we show the same situation depicted in Fig. |5.2| using the Hartree 
approximation. The initial condition on the zero mode is r]{TQ = 0) = -\/3 • 10^''; 
the factor of -\/3 appears due to the different scaling in the zero mode equation. 



( p.89|) , which causes the minimum of the tree level effective potential in the Hartree 



approximation to have a value oi r] = ^/S. Again, the Hubble constant has the 
value h{To) = 0.1. Here, we see again that there is an initial burst of particle 
production as gTi (Fig. p.3| b) grows large. However, the zero mode (Fig. |5.3| a) 

128 



quickly reaches the minimum of the potential and the condition 

-l + r/2(r)/3 + (7S(r) = (5.12) 

is approximately satisfied by forcing the value of gT, quickly to zero. There are 
somewhat fewer particles produced here compared to the large N case, and the 
distribution of particles is more extended. Since the effective mass of the particles is 
nonzero, we expect a matter dominated equation of state (Fig. |5.3| e) for later times. 
The fact that the Hartree approximation does not satisfy Goldstone's theorem 
means that the resulting particles must be massive, explaining why somewhat 
fewer particles are produced. 

Finally, we show the special case in which there is no initial biasing in the field. 



r]{T() = 0) = 0, 'f]{To = 0) = 0, and /i(ro) = 0.1 in Fig. ^]Ja-d. With such an initial 
condition, the Hartree approximation and the large N limit are equivalent. The 
zero mode remains zero for all time, so that the quantity g'^ir) (Fig. |5.4| a) satisfies 
the sum rule (|5.11|) by reaching the value one without decaying for late times. 
Notice that many more particles are produced in this case (Fig. |5.4j b); the growth 
of the particle number for late times is due to the expansion of the universe. The 
particle distribution (Fig. p.4| c) is similar to that of the slow roll case in figure 1. 
The equation of state (Fig. p.4| d) is likewise similar. 

In each of these cases of slow roll dynamics, increasing the Hubble constant 
has the effect of slowing the growth of both rj and gTi. The equation of state 
will be that of a de Sitter universe for a longer period before moving to a matter 
dominated equation of state. Otherwise, the dynamics is much the same as in 
Figs. 
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Figure 5.3: Symmetry broken, slow roll, Hartree, matter dominated evolution of 
(a) the zero mode r){T) vs. r, (b) the quantum fluctuation operator fi'S(T) vs. r, 
(c) the number of particles gN{r) vs. r, (d) the particle distribution gNq{T) vs. 
g at r = 150.7 (dashed line) and r = 396.1 (solid line), and (e) the ratio of the 
pressure and energy density p(r)/£(r) vs. r for the parameter values = — 1, 
r]{ro) = 3V2 . IQ-^, f]{To) = 0, g = lO'^^, h{ro) = 0.1. 
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Figure a 



Figure b 




Figure 5.4: Symmetry broken, no roll, matter dominated evolution of (a) the 
quantum fluctuation operator 5'S(r) vs. r, (b) the number of particles gN{r) vs. r, 
(c) the particle distribution gNqij) vs. q t — 150.1 (dashed line) and r = 397.1 
(solid line), and (d) the ratio of the pressure and energy density p(r)/£(T) vs. r 
for the parameter values — —1, 77(to) = 0, 77(to) — Q,g — 10~^^, /i(to) = 0.1. 
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Case 2: < 0, rj^To) ^ 1. We now examine the case of preheating occurring 
in a chaotic inflationary scenario with a symmetry broken potential. In chaotic 
inflation, the zero mode begins with a value rjlr) ^ 1. During the de Sitter phase, 
h ^ 1, and the field initially evolves classically, dominated by the first order 
derivative term appearing in the zero mode equation [see ( p.89| )]. Eventually, the 
zero mode rolls down the potential, ending the de Sitter phase and beginning the 
preheating phase. We consider the field dynamics in the FRW universe where 
preheating occurs after the end of inflation. We thus take the initial temperature 
to be zero, T = 0. 

Figure ^]5| shows our results for the quantities, rj^r), gT^^r), gN^r), gNg^r), 
and p{T)/e{T) for the evolution in the large N approximation within a radiation 
dominated gravitational background with /i(ro) = 0.1. The initial condition on 
the zero mode is chosen to have the representative value ri{To = 0) = 4 with 
?7(ro = 0) = 0. Initial values of the zero mode much smaller than this will not 
produce significant growth of quantum fluctuations; initial values larger than this 
produces qualitatively similar results, although the resulting number of particles 
will be greater and the time it takes for the zero mode to settle into its asymptotic 
state will be longer. 



We see from Fig. ^.5| a that the zero mode oscillates rapidly, while the am- 
phtude of the oscillation decreases due to the expansion of the universe. This 
oscillation induces particle production through the process of parametric amplifl- 
cation (Fig. p3| c) and causes the fluctuation gT, to grow (Fig. [STSlb). Eventually, 
the zero mode loses enough energy that it is restricted to one of the two minima 
of the tree level effective potential. The subsequent evolution closely follows that 
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Figure 5.5: Symmetry broken, chaotic, large N, radiation dominated evolution of 
(a) the zero mode r]{T) vs. r, (b) the quantum fluctuation operator 5'S(r) vs. r, 
(c) the number of particles gN^r) vs. r, (d) the particle distribution gNq{T) vs. 
q at T = 76.4 (dashed line) and r = 392.8 (solid line), and (e) the ratio of the 
pressure and energy density ]9(r)/e(r) vs. r for the parameter values = — 1, 
r/(ro) = 4, 7)(ro) = 0, ^ = 10-l^ /i(to) = 0.1. 
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of Case 1 above with gT, decaying in time as l/a^(r) ~ 1/r with rj given by the 
sum rule ( |5.11| ). The spectrum (Fig. pl5|d) indicates a single unstable band of par- 



ticle production dominated by the modes g = 1/2 to about q = 3 for late times. 
The structure within this band becomes more complex with time and shifts some- 
what toward lower momentum modes. Such a shift has already been observed in 
Minkowski spacetimes |^ . Figure |5.5| e shows the equation of state which we see 



to be somewhere between the relations for matter and radiation for times out as 
far as r = 400, but slowly moving to a matter equation of state. Since matter red- 
shifts as l/a^(r) while radiation redshifts as l/a^(r), the equation of state should 
eventually become matter dominated. Given the equation of state indicated by 
Fig. |5.5| e, we estimate that this occurs for times of order r = 10^. The reason 
the equation of state in this case differs from that of cold matter as was seen in 
Figs. |5.2|-|5.4| is that the particle distribution produced by parametric amplification 



is concentrated at higher momenta, g ~ 1. 

Figure |5.6| shows the corresponding case with a matter dominated background. 



The results are qualitatively very similar to those described for Fig. ^.5| above. Due 
to the faster expansion, the zero mode (Fig. p.6| a) finds one of the two wells more 
quickly and slightly less particles are produced. For late times, the fluctuation 
gT, (Fig. ^.6| b) decays as l/a^(r) oc l/r^^^. Again we see an equation of state 
(Fig. f).6^ ) which evolves from a state between that of pure radiation or matter 
toward one of cold matter. 

The Hartree case is depicted in Fig. |5.7| for a matter dominated universe, with 
the initial condition on the zero mode 'r](TQ = 0) = A^/S. Again, the evolution 
begins in much the same manner as in the large N approximation with oscillation 
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Figure 5.6: Symmetry broken, chaotic, large N, matter dominated evolution of (a) 
the zero mode r){T) vs. r, (b) the quantum fluctuation operator fi'S(r) vs. r, (c) 
the number of particles gN^r) vs. r, (d) the particle distribution gNq{T) vs. q at 
r = 50.8 (dashed line) and r = 399.4 (solid line), and (e) the ratio of the pressure 
and energy density p{T)/e{T) vs. r for the parameter values = —1, ?7(ro) = 4, 
v{ro) = 0,g = 10-^^ h{To) = 0.1. 
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of the zero mode (Fig. |5.7| a), which eventually settles into one of the two minima 
of the effective potential. Whereas in the large approximation, the zero mode 
approaches the minimum asymptotically [as given by and our late time 

analysis below], in the Hartree approximation we see that the zero mode finds 
the minimum quickly and proceeds to oscillate about that value. The two point 
correlator (Fig. |5.71 b) quickly evolves toward zero without growing large. Particle 
production in the Hartree approximation (Figs. |5.71 c-d) is again seen to be ineffi- 



cient compared to that of the large N case above. Figure ^.7| 8 again shows that 
the equation of state is matter dominated for all but the earliest times. 

A larger Hubble constant prevents significant particle production unless the 
initial amplitude of the zero mode is likewise increased such that the relation 
?7(ro) ^ h{To) is satisfied. For very large amplitude ^^(to) ^ 1, to the extent that 
the mass term can be neglected and while the quantum fluctuation term has not 
grown to be large, the equations of motion ( |3.89D and ( |3.9CID are scale invariant 
with the scaling r] — > fn], h — » ^h, r ^ r/fi, and q — > fiq, where is an arbitrary 
scale. 

For completeness, we show the case of the evolution with initial values of the 



Hubble constant given by h(To) = 5 and /^(to) = 2 respectively in Figs. |5.8| and 



using radiation dominated expansion. Here, we have used the large N approxima- 
tion and have made an appropriate increase in the initial value of the zero mode 
such that the fluctuations (Figs. |5.^ and |5.9| b) grow significantly [we have chosen 
vi'^o) = 40 in Fig. |5.8| and tjItq) = 16 in Fig. ^.9|| . While the dynamics looks 
much hke that of Figs. |575|-|5r7| above, we point out that the particle distribution 
(Figs. p78| d and |0|d) is extended to higher values with the result being that the 
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Figure 5.7: Symmetry broken, chaotic, Hartree, matter dominated evolution of 
(a) the zero mode r]{T) vs. r, (b) the quantum fluctuation operator 5'S(r) vs. r, 
(c) the number of particles gN{r) vs. r, (d) the particle distribution gNq{T) vs. 
q at T = 151.3 (dashed line) and r = 397.0 (solid line), and (e) the ratio of the 
pressure and energy density p{T)/e{T) vs. r for the parameter values = — 1, 
r/(ro) = 4 ■ 31/2, f,{To) = 0, g = IQ-'^ h{ro) = 0.1. 
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equation of state (Figs. |5.8| e and ^.Ql e) is weighted more toward that of radiation. 

Case 3: > 0, //(tq) ^ 1. The final case we examine is that of a simple 
chaotic scenario with a positive mass term in the Lagrangian. Again, preheating 
can begin only after the inflationary phase of exponential expansion; this allows 
us to take a zero temperature initial state for the FRW stage. 



Figure |5.10| shows this situation in the large N approximation for a matter 
dominated cosmology. The zero mode, f]{T), oscillates in time while decaying in 
amplitude from its initial value of ri{TQ = 0) = 5, r]{TQ = 0) = (Fig. |5.1(J| a), while 
the quantum fluctuation, gT,, grows rapidly for early times due to parametric 
resonance (Figs. |5.1CI| b). We choose here an initial condition on the zero mode 



which differs from that of Figs. |5.5| - |5.6| above since there is no significant growth 
of quantum fluctuations for smaller initial values. From Fig. 5.10| d, we see that 
there exists a single unstable band at values of roughly g = 1 to g = 3, although 
careful examination reveals that the unstable band extends all the way to g = 0. 
The equation of state is depicted by the quantity p(r)/£:(r) in Fig. |5.1CI| e. As 
expected in this massive theory, the equation of state is matter dominated. 



The final case is the Hartree approximation, shown in Fig. Here, paramet- 
ric amplification is entirely inefficient when expansion of the universe is included 
and we require an initial condition on the zero mode of r]{TQ = 0) = 12^/3 to 
provide even meager growth of quantum fluctuations. We have used a matter 
dominated gravitational background with /i(ro) = 0.1. We see that while the 
zero mode oscillates (Fig. |5.11| a), there is little growth in quantum fluctuations 
(Fig. |5.11[ b) and few particles produced (Fig. |5.11| c). Examining the particle dis- 
tribution (Fig. |5.11| d), it is found that the bulk of these particles is produced within 
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Figure 5.8: Symmetry broken, chaotic, large N, radiation dominated evolution of 
(a) the zero mode r]{T) vs. r, (b) the quantum fluctuation operator 5'S(r) vs. r, 
(c) the number of particles gN{r) vs. r, (d) the particle distribution gNq{T) vs. 
q at T = 118.9 (dashed line) and r = 394.7 (solid line), and (e) the ratio of the 
pressure and energy density p(r)/e(r) vs. r for the parameter values = — 1, 
r/(ro) = 40, f]{To) = 0,9 = IQ-'^ h{ro) = 5.0. 
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Figure 5.9: Symmetry broken, chaotic, large N, radiation dominated evolution of 
(a) the zero mode r){T) vs. r, (b) the quantum fluctuation operator fi'S(T) vs. r, 
(c) the number of particles gN{r) vs. r, (d) the particle distribution gNq{T) vs. 
g at r = 55.1 (dashed line) and r = 194.2 (solid line), and (e) the ratio of the 
pressure and energy density p(r)/e(r) vs. r for the parameter values = — 1, 
r]{ro) = 16, 7)(ro) = 0, ^ = IQ-^^, h{To) = 2.0. 
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Figure 5.10: Symmetry unbroken, chaotic, large A^, matter dominated evolution 
of (a) the zero mode tjIt) vs. r, (b) the quantum fluctuation operator g'S(T) vs. 
T, (c) the number of particles gN{r) vs. r, (d) the particle distribution gNqij) 
vs. g at r = 77.4 (dashed line) and r = 399.7 (solid line), and (c) the ratio of the 
pressure and energy density ]9(r)/e(r) vs. r for the parameter values w? = +1, 
riiro) = 5, 7)(ro) = 0, g ^ lO'^^ h{To) = 0.1. 
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a single resonance band extending from 15tog~ 16. This resonance develops 
at early time during the large amplitude oscillation of the zero mode. These results 
are explained by a simple resonance band analysis described below. 

At first glance, it is not entirely clear why there are so many more particles 



produced in the large N case of Fig. 5.10 than in the Hartree case of Fig. 5.11 



Since in the present case the Hubble time is long compared to the oscillation time 
scale of the zero mode, h <^ 1, we would expect a forbidden band for early times at 



the location given approximately by the Minkowski results provided in Ref. [Q. 
In fact, we find this to be the case. However, we know from previous studies that 
in Minkowski space a similar number of particles is produced in both the Hartree 



and large case [l54| , pS 



The solution to this problem is inherent in the band structure of the two cases 
when combined with an understanding of the dynamics in an expanding spacetime. 
First, we note that, for early times when gH <^ 1, the zero mode is well fit by 
the function ?7(r) = ?7o/(r)/a(r) where /(r) is an oscillatory function taking on 
values from —1 to 1. This is clearly seen from the envelope function rio/a{T) 
shown in Fig. 5.11| a (recall that gJ] <^ 1 during the entire evolution in this case). 



Second, the momentum that appears in the equations for the modes (|3.90|) is the 
physical momentum q/a{T). We therefore write the approximate expressions for 
the locations of the forbidden bands in FRW by using the Minkowski results of 
5^ with the substitutions t/q — > '-)7]Q/a?{T) (where the factor of 7 accounts for the 



difference in the definition of the non-linear coupling between this study and 
and — > /a?{T). 
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Figure 5.11: Symmetry unbroken, chaotic, Hartree, matter dominated evolution 
of (a) the zero mode 'r]{T) vs. r, (b) the quantum fluctuation operator g'S(T) vs. 
T, (c) the number of particles gN{r) vs. r, (d) the particle distribution gNqij) 
vs. g at r = 50.5 (dashed line) and r = 391.2 (solid line), and (e) the ratio of the 
pressure and energy density p{t)/s{t) vs. r for the parameter values = +1, 
riiro) = 12 • 3V2, fi{To) =0,g^ IQ-^^, /^(ro) = 0.1. 
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Making these substitutions, we find for tlie location in comoving momentum q 
of the forbidden band in the large N (Fig. p.lO| ) and Hartree (Fig. cases: 

I 
2 



< < — , (large N) (5.13) 



+ + 4 + 1 I . (Hartree) (5.14) 



The important feature to notice is that while the location of the unstable band (to 
a first approximation) in the case of the continuous 0{N) theory is the same as in 
Minkowski and does not change in time, the location of the band is time dependent 
in the discrete theory described by the non-perturbative Hartree approximation. 
While rio/a{r) ^ 1, the Hartree relation reduces to 

|<,^<|. (5.15) 

This is the same as the Minkowski result for large amplitude, and one finds that 
this expression accurately predicts the location of the resonance band of Fig. |5.11| d. 
However, with time the band shifts its location toward higher values of comoving 
momentum as given by ( |5.14D , cutting off particle production in that initial band. 



There is continuing particle production for higher modes, but since the Floquet 
index is decreased due to the reduced amplitude of the zero mode, since there is no 
enhancement of production of particles in these modes (as these modes begin with 
at most of order 1 particles), and because the band continues to shift to higher 
momenta while becoming smaller in width, this particle production never becomes 
significant. 

We emphasize that this result is not an artifact of the approximations used but 
rather reflects an important difference between the behaviors of the theory with 

144 



discrete symmetry described by the Hartree approximation and the theory with 
continuous symmetry described by the large hmit. In the case of a continuous 
symmetry, the dynamics must satisfy the appropriate Ward identities correspond- 
ing to that symmetry. The Hartree approximation does not respect these identities 
and is therefore only appropriate for describing the discrete symmetry. Likewise, 
the discrete case is not well described by the large N approximation. Mathemati- 
cally, these differences show up in the form of the differing band structure described 
above, and we conclude that preheating is inefficient in simple models in which the 
relevant field obeys an unbroken discrete symmetry. 

As in the symmetry broken case of Figs. ^.5| - 0.7|, the equations of motion for 



large amplitude and relatively early times are approximately scale invariant. In 
Figs. 5.12| - 5.13 we show the case of the large N evolution in a radiation dominated 



universe with initial Hubble constants of /i(ro) =5 and /?.(ro) =2 respectively with 
appropriately scaled initial values of the zero mode of i]{tq) = 40 and //(tq) = 16. 
Again, the qualitative dynamics remains largely unchanged from the case of a 
smaller Hubble constant. 

5.2.3 Late Time Behavior 

We see clearly from the numerical evolution that in the case of a symmetry broken 



potential, the late time large solutions obey the sum rule ( |4.7| ). Because of this 
sum rule, we can write down the analytical expressions for the late time behavior of 
the fluctuations and the zero mode. Using ([4.7|), the mode equation ( p.90| ) becomes 



a(t) d P 



dt^ a{t)dt a?{t) 



fk{t) = 0. (5.16) 
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Figure 5.12: Symmetry unbroken, chaotic, large N, radiation dominated evolution 
of (a) the zero mode r7(r) vs. r, (b) the quantum fluctuation operator g'S(T) vs. 
T, (c) the number of particles gNi^r) vs. r, (d) the particle distribution gNq{T) vs. 
g at r = 117.3 (dashed line) and r = 393.6 (solid line), and (e) the ratio of the 
pressure and energy density p(r)/£(r) vs. r for the parameter values = +1, 
r]{ro) = 40, 7)(ro) = 0, ^ = IQ-^^, h{To) = 5.0. 
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Figure 5.13: Symmetry unbroken, chaotic, large N, radiation dominated evolution 
of (a) the zero mode tjIt) vs. r, (b) the quantum fluctuation operator g'S(T) vs. 
r, (c) the number of particles gN^r) vs. r, (d) the particle distribution gNq{T) vs. 
g at r = 102.1 (dashed line) and r = 251.6 (solid line), and (e) the ratio of the 
pressure and energy density p(r)/e(r) vs. r for the parameter values = +1, 
r/(ro) = 16, f]{To) = 0, g = IQ-'^ h{ro) = 2.0. 
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This equation can be solved exactly if we assume a power law dependence for the 
scale factor a(t) = (t/to)"- The solution is 

+ 4 t^'-''^^" , (5.17) 

where and Y^, are Bessel and Neumann functions respectively, and the constants 
Cfc and dk carry dependence on the initial conditions and on the dynamics up to 
the point at which the sum rule is satisfied. 

These functions have several important properties. In particular, in radiation 
or matter dominated universes, n < 1, and for values of wavenumber satisfying 
k ^ t^^^^^'^tQ, the mode functions decay in time as l/a(t) ~ t"". Since the sum 
rule applies for late times, r — tq ^ 1, we see that all values of k except a very small 
band about k = redshift as l/a{t). The k = mode, however, remains constant in 
time, explaining the support evidenced in the numerical results for values of small 
k (see Figs. and p^ ). These results mean that the quantum fluctuation has a 
late time dependence of {ip^it))^ ~ l/a^{t). The late time dependence of the zero 
mode is given by this expression combined with the sum rule (|4.7| ). These results 
are accurately reproduced by our numerical analysis. Note that qualitatively this 
late time dependence is independent of the choice of initial conditions for the zero 
mode, except that there is no growth of modes near /c = in the case in which 
particles are produced via parametric amplification (Figs. |5.5| and ^.6|) . 

For the radiation n = ^ and matter dominated ^ = | universes, Eq. ( |5.17] ) 
reduces to elementary functions: 

a{t)Mt) = c, e'<''''^' + d, e-'<'''^' (RD) , 
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«(^) fk{t) = Ck € 




1 + 



I 



+ 4 e 




I 



(MD) . 



(5.18) 



It is also of interest to examine the n > 1 case. Here, the modes of interest 
satisfy the condition k <C t"~^/to for late times. These modes are constant in time 
and one sees that the modes are frozen. In the case of a de Sitter universe, we can 
formally take the limit n ^ oo and wc sec that all modes become frozen at late 
times. This case was studied in detail in Ch. 4. 

5.3 Conclusions 

We have shown that there can be significant particle production through quantum 
fluctuations after inflation. However, this production is somewhat sensitive to the 
expansion of the universe. From our analysis of the equation of state, we see 
that the late time dynamics is given by a matter dominated cosmology. We have 
also shown that the quantum fluctuations of the inflaton decay for late times as 
l/a^(i), while in the case of a symmetry broken inflationary model, the inflaton 
fleld moves to the minimum of its tree level potential. The exception to this 
behavior is the case when the inflaton begins exactly at the unstable extremum of 
its potential for which the fluctuations grow out to the minimum of the potential 
and do not decay. Initial production of particles due to parametric amplification 
is significantly greater in chaotic scenarios with symmetry broken potentials than 
in the corresponding theories with positive mass terms in the Lagrangian, given 
similar initial conditions on the zero mode of the infiaton. 

Since there are a number of articles in the literature treating the problem of 
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preheating, it is useful to review the unique features of the present work. First, 
we have treated the problem dynamically, without using the effective potential 
(an equilibrium construct) to determine the evolution. Second, we have provided 
consistent non-perturbative calculations of the evolution to bring out some of the 
most relevant aspects of the late time behavior. In particular, we found that 
the quantum backreaction naturally inhibits catastrophic growth of fluctuations 
and provides a smooth transition to the late time regime in which the quantum 
fluctuations decay as the zero mode approaches its asymptotic state. Third, the 
dynamics studied obeys the constraint of covariant conservation of the energy 
momentum tensor. 

As we have seen, significant particle production may occur due to parametric 
amplification after inflation. However, this production depends strongly on the 
expansion rate and is effectively cut off for values of the expansion rate h{to) > 1. 
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Chapter 6 
Conclusions 



Presently, inflation is the most promising paradigm for solving the horizon, flatness, 
and relic problems, as well as the fluctuation problem. However, it remains a 
paradigm in search of foundation with literally thousands of models, none of which 
are particularly well motivated by an overall theory of particle physics. As new 
observational data on cosmological parameters become available, we should begin 
to see the elimination of large classes of these models, and there is certainly the 
potential that a small class of models will gain favor. The ideal situation, of course, 
would be for such clues as to the nature of inflation to provide guidance as to the 
nature of particle physics in general at energies approaching the Grand Unifled 
scale. 

This is an ambitious goal and will not only require unprecedented observational 
tools such as the MAP and PLANCK satellites, but also complete utilization of 
all available theoretical tools. Observations may provide values of cosmological 
parameters to an accuracy of one percent within the next decade, and it will 
be necessary to have confidence in theoretical predictions of these parameters for 
various cosmological models to a similar accuracy. It therefore becomes necessary 
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to question the methods that have been utihzed during the past 20 years as these 
methods were required to provide only order of magnitude predictions. 

Hence, we have undertaken this study of the non-equihbrium aspects of in- 
flationary theory, and we have discovered a number of important properties of 
inflationary dynamics. Of particular note is the construction of a classical back- 
ground field from the quantum inflaton in new inflation models, which allows us 
to make confident predictions of the primordial density perturbations. Also of 
primary importance is the discovery that explosive particle production may not 
be significant if the expansion rate immediately after inflation is of the order or 
larger than the inflaton mass scale. In addition, we have shown how the symmetry 
properties of the underlying inflationary model may have a significant influence on 
the importance of explosive particle production in an expanding spacetime. 

While we have used the powerful tools of non-equilibrium field theory to study 
one of the most exciting areas in cosmology, we note that there are a number of 
other problems to which these techniques are particularly well suited. Baryogenesis 
is a striking example, since it cannot take place at all in equilibrium. A thorough 
out of equilibrium study is bound to provide new insight into this process, which 
results in the production of the bulk of the matter with which we are familiar. 

Another process taken from cosmology is topological defect formation. Topo- 
logical defects represent the only known competition to inflation as the source of 
primordial density perturbations. And while they are presently out of favor as 
candidates for seeding structure, they will remain important because they are al- 
most certain to have formed during the evolution of the universe. They may be a 
source of dark matter, or the necessity of mechanisms which reduce their number 
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and importance may provide important clues about the history of our universe. 

The potential uses of non-equilibrium field theory outside of cosmology are 
just as numerous, with many applications in condensed matter physics as well as 
particle physics, most of which remain unexplored. 

This is an exciting time for cosmology and physics. We are now just beginning 
to discover the answers to some of the most profound questions that can be asked: 

• Is the universe open or flat? 

• What kind of physics is there beyond the Standard Model of Particles? 

• What is the nature of the dark matter? 

• Is the universe accelerating? 

• What would it mean to have a significant cosmological constant? 

• How old is the universe? 

• What does this all mean for the physics of the early universe? 
There is only one thing for certain: there will be surprises! 
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